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ABSTRACT 

We use the high spatial and spectral resolution of the PAWS CO (1-0) survey of the 
inner 9 kpc of the iconic spiral galaxy M51 to examine the effect of gas streaming mo- 
tions on the star-forming properties of individual GMCs. We compare our view of gas 
flows in M51 - which arise due to departures from axi-symmetry in the gravitational 
potential (i.e. the nuclear bar and spiral arms) - with the global pattern of star forma- 
tion as traced by Ha and 24/xm emission. We find that the dynamical environment of 
GMCs strongly affects their ability to form stars, in the sense that GMCs situated in 
regions with large streaming motions can be stabilized, while similarly massive GMCs 
in regions without streaming go on to efficiently form stars. We argue that this is the 
result of reduced surface pressure felt by clouds embedded in an ambient medium un- 
dergoing large streaming motions, which prevents collapse. Indeed, the variation in gas 
depletion time expected based on the observed streaming motions throughout the disk of 
M51 quantitatively agrees with the variation in observed gas depletion time scale. The 
example of M51 shows that streaming motions, triggered by gravitational instabilities in 
the form of bars and spiral arms, can alter the star formation law; this can explain the 
variation in gas depletion time among galaxies with different masses and morphologies. 
In particular, we can explain the long gas depletion times in spiral galaxies compared 
to dwarf galaxies and starbursts. We suggest that adding a dynamical pressure term to 
the canonical free-fall time produces a single star formation law that can be applied to 
all star-forming regions and galaxies, across cosmic time. 



1. INTRODUCTION 

Our understanding of how galaxies build up 
their mass - fundamental to our picture of galaxy 
evolution - relies heavily on the relation between 
newborn stars and the dense molecular material 
out of which they form. At p resent, the em- 
pirical Keimicu tt-Schmidt (KS; I Schmidt I 119591 : 
iKennicuTa 119891 ) relation linking the star forma- 
tion rate surface density ^sfr with a power- 
law dependence on gas surface density E gas is in 
wide use by numerical simulations that inform our 
view of galaxy evolution. This relation has been 
demonstrated to hold over several orders of mag- 
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nitude in S gas and across a range of galaxy types. 
The first homogeneous spatiall y-resolved observa- 
tions in nearby galaxies (e.g. iBigiel et al.l [2008: 
Leroy et al. 2008,2009,2012) suggest that gas is 
consumed in star formation at a fixed efficiency, 
leading to an apparently 'universal', but long, gas 
depletion time of ~2.5 Gyr for galaxies. The rela- 
tion was further shown to be determined predom- 
inant ly by the molecular phase of the n eutral gas 
(|Schruba et al.ll2lrilllBTgiel et alJ[2008]) . 

At the same time, there is mounting evidence 
that star formation may occur in two modes: 
the 'normal' mode typical of the majority of 
star formation in local disks, contrasted by a 
'starburst' mode, wherein gas is consumed much 
more efficiently and quickly (i.e. iKennicutva 
1998). Bimodality in the star format ion (SF) re- 
lation is manifest both local ly (i.e. iLerov et al.l 
120131 : iSaintonge et~all 120121) and at hi g her z 
(i.e. iDaddi et al.| 120101 : lGenzel"et~aT1 [Mot 



Garcia-Burill o et aLll2012t ). But, to date, there is 
no single 'universal' star formation relation than 
can smo othly link the two regimes, across a ll times 
(but see lKrumholz. Dekel fc McKee ||201~2| ). 

How the star-forming gas is organized may in- 
fluence the depletion time (jCombes et al.l [2012; 
Saintong e et al.ll2012l and references therein). An 
interpretation of the long depletion time in 'nor- 
mal' star forming galaxies could be that it re- 



fleets only the (low) fraction of star-forming 
clouds rather than the fraction of a cloud that 
forms stars. Variations in the properties of Gi- 
ant Molecular Clouds (GMCs), the actual seeds 
of star formation, themselves could profoundly 
impact their star formation capabilities as ev- 
idenced by the large scatter of observed star 
formation rates (SFRs) in lo cal Galactic GMCs 
(iMoonev fe Solomon! 119881: iKennicutt fe Evansl 
I2012t iHeiderman et all I2010T ). Yet, our under- 
standing of the molecular gas phase and its organi- 
zation into self-gravitati ng entities, i.e. GMCs, i s 
limited (for a review, see McKce & Ostriker 2007). 

Observations of molecular gas at high reso- 
lution in nearby spiral galaxies are critical to 
test the emerging picture of how galactic envi- 
ronment influences the organization, properties 
and evolutionary pathways of a galaxy's GMC 
population, and its global patterns of star for- 
mation. In this paper we continue our study 
of the interacting Whirlpoo l galaxy M51 (D=7.6 
Mpc; ICiardullo et al.l [20021 ). which is more rep- 
resentative of a typical star forming galaxy than 
the well-studied lower mass galaxies in the Lo- 
cal Group. The unparalleled high resolution of 
the PAWS (PdBI Arc second Whirlpool Survey; 
Schinncr er et al.l I2013D data, combined with ex- 
ceptional multi-wavelength coverage, makes it an 
ideal target for examining the influence of bar and 
spiral instabilities - and more generally, dynami- 
cal galactic environment - on the organization of 
the ISM and global star formation patterns. 

Our detailed analysis of the molecular gas in 
M51- including its GMC population - so far sug- 
gests that seve ral mechanisms control its structure 
and evolution ([Colombo et al.ll2~013a[ ). Deviations 
in the gravitational potential from axisymmetry 
(i.e. the nuclear bar and two-armed spiral) induce 
shear, shocks, and streaming motions in the gas 
flow, leading to changes in the local gas surface 
density. This subsequently alters the pattern of 
radiation from recent star formation, which oc- 
curs on the convex side of the spiral arm, and 
therefore also the mechanical energy input from 
newly formed OB stars, stellar winds and super- 
novae. The galactic environment is thus very dy- 
namic in both senses, and it is not clear under 
which conditions the developed static models ap- 
ply (e.g. Krumholz et al. 2009). 

Gas kinematics on the scales of individual 
GMCs are key for linking on-going star forma- 
tion to the organization and radial transport of 
gas by the stellar spiral arms. We expect the dy- 
namical environment to strongly impact the orga- 
nization of the ISM in a galaxy like M51, which 
hosts strong, well-organized spiral structure and is 
presently undergoing interaction with a compan- 
ion. Depending on location in the disk, gas will be 
driven radially in- or outward and may occasion- 
ally be hal ted along its path , build ing up rings (see 
review by iButa fc Combesl ([19961 )). Whether gas 
is in motion or stationary will naturally influence 
the build-up of the star-forming gas reservoir, the 
accumulation of mass in to clouds, and the even- 



tual collapse to form stars. 

To quantify the influence of gas kinematics on 
cloud scales we undertake a study in two parts. 
We begin by evaluating the close link between 
the gas depletion time Tdep='Egas/^SFn-J and dy- 
namical environment as defined by the radial pro- 
file of present-day gravitational torques f ^3.1[) . We 
then establish the magnitude and size of gas mo- 
tions driven in response to the spiral torquing 
f£ J3.2p and interpret the behavior in Td ep in terms 
of gas kinematics in § 14.11 We suggest that dy- 
namical pressure plays a role in stabilizing GMCs. 

Second, we develop a picture that links gas 
streaming motions to the gas depletion time using 
the Bernoulli principle, which equates an increase 
in gas velocity with reduced external cloud pres- 
sure (§ 14. 2p . We use the radial profiles of stream- 
ing motions to assess the role of shear on the 
stability of gas and discuss the evidence against 
shock- and stellar feedback-driven turbulent sup- 
port in § 15.11 We conclude with a discussion of 
the broader implications of this result, from the 
angular offset between gas and young star tracers 
to the scatter in the KS star formation relation (§ 
15. 2[) . In Q5.'3\ we present a new form for the star 
formation relation between ^sfr and E gas that 
smoothly links star formation at low-z and high-z 
by parameterizing the environmental dependence 
of clouds in terms of non-circular gas streaming 
motions. In particular, we predict a link between 
gas depletion ti me and gas fraction, as recently ob- 
served locally (iSaintonge et alj|2012l) and in evo- 
lution from z=2 (Comb es et al.ll2012j ). 

2. DATA 

The analysis in this paper is based on CO (1-0) 
data obtained by the Plateau de Bure Arcsec- 
ond Whirlpool Survey (PAWS, iSchinnerer et al.l 
2013). PAWS observations cover a roughly 
270" x 170" field of view in the inner disk of M51 
and include both interferometric data from the 
Plateau dc Bure Interferometer (PdBI) and single- 
dish data from the IRAM 30m telescope. The 
effective angular resolution in the final combined 
data cube is 1. "16x0. "97, or a spatial resolution 
of ^40 pc at our assumed distance of M51. The 
channel separation is 5 km s _1 . The PAWS ob- 
serving strategy, data reduction and combination 
procedures , and flux calibration are described by 
IPetvet al.l(l20Tl . 

Here we use maps of the zeroth, first and sec- 
on d moments of the PAWS data cube as presented 
bv IPetvet all (|2013[ ). From the PAWS map of in- 
tegrated CO(1-0) intensities (moment zero), we 
estimate the molecular gas surface density, ap- 
plying a constant conversion facto r Xgo= 2x 10 20 



cm 2 (K km s 1 )~ 

10 In the extragalactic context, the gas depletion time 
is equated with the inverse of the star formation efficiency 
SFE=Sgi?fl/S 9as , although the former quantity is not a 
traditional efficiency, but an efficiency per free fall time. 
This is in contrast with the star formation efficiency e re- 
ferred to in galactic studies, which represents the fraction 
of gas converted in to stars per star formation event. 
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Fig. 1. — Radial profiles of the azimuthally-averaged 
torque (in units of the absolute magnitude of the aver- 
age torque across the PAWS field of view <T>|; black 
and white) and molecular gas depletion time AT,j e p=Td e p- 
<'Tde.p> (where <Td ep > is the average t^j, across the 
PAWS field of view; blue) measured in 2.4" radial bins. 
Each crossing from negative to positive torque corresponds 
to the location of the corotation radius (CR) of the struc- 
ture: inside CR material is driven radially inward and out- 
side material moves outward. 



20133) an( A including a factor of 1.36 to account 
for the presence of Helium. From the line-of-sight 
velocity field (first moment) we derive estimates 
for the circular velocity and the non-circular mo- 
tions. The velocity dispersion map (moment two) 
supplies our measure of the turbulent gas mo- 
tions after removal of the contribution from un- 
resolved bulk motions estimated from the line-of- 
sight vel ocity field (adop ting the formalism devel- 
oped bylPetv et aLll2013|) . We refer the reader to 
iPetv et al.l (|2013t ) for a study of the resolution- 
dependence of observables such as the maximum 
CO brightness (i.e. gas surface density) and veloc- 
ity second moment (gas velocity dispersion) mea- 
sured from the PAWS data cube. 

The PAWS data are complemented by an exten- 
sive multi- wavelength dat a set, which has been as - 
sembled and described in lSchinnerer et a l. (2013). 
Here we make use of the S 4 G IRAC 3.6 /im image 
(|Sheth et al.ll20lll) . cor rected for the cont ribution 
of non-stellar emission (|Meidt et al.ll2012h . to de- 
fine the stellar mass distribution. As described in 
Appendix IB") we use this stellar mass map to de- 
rive the stellar potential and the profile of present 
day torques. We also use the Ha image (corrected 
for the stellar contin uum using the procedure of 
iGutierrez et al.l 1201 ID an d MIPS 24 urn data ac- 
quired as part of SINGS (jKennicutt et al. I [20031 
to trace the star formation in M51. The MIPS 
24 fim image has been processed with the HiRes 
algorith m, giving a lim iting angular resolution of 
2.4" l|Dumas e t al. 201lJ). We combine the two im- 
ages to trace the total ( obscured and unob scured) 
star formation following iLerov et al.1 (120091) usin g 
the empirical calibration of ICalzetti et al.l (|2007). 



3. GAS FLOW AND GLOBAL STAR 
FORMATION PATTERNS IN M51 

M51 is a favorite test-bed for spiral arm den- 
sity wave theories. Studies of the spiral morphol- 
ogy and kinematics show evidence for the offset 
alignment of the gaseous, young and old stellar 
tracers predicted by theory, as well as strong non- 
circular gas st reaming mo tions ([Vogel et al.lll 99" 
Garcia-Buril lo et al.lll9 93: Rand 19 93tlAalto et al 
199 9!: ISchuster et al.1 I2007t IShettv et al.l I2007t 
Schinncrcr et al. 2013). In this section, we ex- 
amine gas motions in M51 using two independent 
approaches. We calculate the present-day torques 
on the molecular gas due to the stellar compo- 
nent. We then decompose the CO velocity field, 
and show that the gas motions implied by our 
kinematic analysis are in excellent agreement with 
the response that we would expect based on our 
torque map. Finally, we compare our picture of 
gas flows with the global pattern of star forma- 
tion. 

3.1. Present-day Torques 

As described in detail in Appendix [Bj our 
2D map of the ol d stellar light in M51 
(|Schinnerer et al.|[2013l) gives us a snap-shot of the 
torques exerted by the non-axisymmetric struc- 
ture (i.e. nuclear bar, spiral arms) present in 
the density and hence gravitational potential of 
the system. These torques drive radially in- and 
outward gas motions, depending on the sign of 
torque. Crossings from negative to positive in 
the azimuthally averaged torque profile (Figure 
[1]) mark the location of corotation, which is de- 
fined as the radius where the angular speeds of 
the disk and non-axisymmetric structure are the 
same. Inside (outside) corotation, negative (posi- 
tive) torques drive gas radially inward (outward). 
Crossings from positive to negative torque coin- 
cide with the switch in predominance to a unique 
pattern. 

This series of zero-crossings defines a set of dy- 
namical environments (see Appendix |B|) : the zone 
dominated by the nuclear bar (i?<20"), the molec- 
ular ring (i?-30"), the inner spiral (30"< J R<85") 
and the outer spiral (i?>85"). The inner spiral 
can be further split in to two zones, inside and 
outside the main spiral's corotation radius (CR) 
at R=60". In what follows, we use this series of 
dynamical environments - and the torques within 
them - as the basis for interpreting the observed 
pattern of non-circular gas streaming motions in 
the inner disk of M51. 

3.2. Radial and azimuthal streaming motions 

If our profile of the present-day torques in M51's 
inner disk is correct, then we expect the signa- 
tures of the gas motions that occur in response 
to the torques to be evident in the velocity field 
derived from the CO data. In Appendix \C\ we 
use our mass-based model of the circular motion 
in M51 (Appendix IA1 to isolate the contribution 
of non-circular motion to the line-of-sight velocity 
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Fig. 2. — (Top) CO intensity-weighted radial (black) and 
azimuthal (gray) non-circular streaming velocities mea- 
sured in 1.5" radial bins, reconstructed as described in the 
text. The intensity- weighting focuses our measurements 
on the CO-bright spiral arms. Thin lines represent the 
rms dispersion in solutions comprised of ±5° variation in 
the assumed spiral arm pitch angle i p added in quadrature 
with the uncertainty introduced by our adopted rotation 
curve (see text). (Bottom) Magnitude of the non-circular 

streaming velocities \v$\ = Jvlj, + v\ (black), derived from 

the non-circular motions above. Thin lines show the uncer- 
tainty measured from the errors above. The gas velocity 
dispersion profile cr-calculated as the azimuthal average of 
the line-of-sig ht dispersion corr ected for unresolved bulk 
motions (i.e. Pcty ct al. 2013)— is shown for comparison 
(dotted). 



field of the CO emission. Then we decompose this 
non-circular motion into its radial and azimuthal 
components v r and v^. 

Figure [2] shows radial profiles of both recon- 
structed velocity components, which are calcu- 
lated using CO intensity as a weighting factor in 
binned azimuthal averages to reveal the velocities 
characteristic of the (CO bright) spiral arms. Pro- 
files with inverse-intensity weighting, which high- 
light the inter-arm streaming velocities, are quali- 
tatively similar but show velocities ~ 5-10 km s _1 
lower than the profiles in Figure [2] (see Appendix 

Our reconstructed radial velocities present us 
with a picture of radial gas flow that is quali- 
tatively similar to that implied by the present- 
day torques. We once again retrieve the pat- 
tern of radial inflow and (weak) outflow that 
mark the locations of at least two distinct CR, 
at nearly the same radii as found with the torque 
map (Figure [5]). Radial velocities go to zero 
and switch sign near 25" and again between the 
range 60"<i?<80" as predicted by the torques. 



Meanwhile, the azimuthal velocities exhibit sim- 
ilar zero-crossings as well as a pattern of posi- 
tive and negative values that reflect a complex 
response to the spiral arm; depending on location 
relative to the spiral, gas at some points is station- 
ary in azimuth or even flows opposite to the direc- 
tion of rotation. There is some dependence on the 
adopted rotation curve (and hence stellar M/L), 
which shifts the azimuthal velocity horizontally in 
Figure[21 But, as described below, this choice does 
not greatly alter the measured magnitude of the 
streaming motions. 

Almost independent of the particular rotation 
curve (or adopted stellar M/L; see Appendix ICl) . 
the regions of negative torque identified in Fig- 
ure [1] map to regions of the highest non-circular 
velocities in the arm. The magnitude of stream- 
ing motions \v$\ = \/v 2 . + v 2 , is shown in the bot- 
tom panel of Figure [2] for each of our three pos- 
sible rotation curves, i.e. assuming the nominal 
M/L=0.45±0.15 (chosen as described in Appendix 
ICJ) . Although we do not expect negative torques 
to coincide with the largest streaming motions as 
a rule, the link between v r and vs is expected, 
according to the predicted relation between radial 
and azimuthal velocities in response to a spiral 
poten tial perturbation (i.e. IBinnev fc Tremaind 
Il987t ): for an m-armed spiral rotating at pat- 
tern speed 51p in a galaxy with a flat rotation 
curve (zero background shear) v ( p=v r 2n/m(fl — 
H, p ) (with il the disk angular velocity), and so the 
total magnitude of streaming motions will follow 
v r , even at CR since v r (x m(VL — fi p ). 

3.3. A Link Between Gas Flows and Star 
Formation Efficiency in M51 's inner disk 

A number of previous studies have shown that 
the pattern of M51's star for mation exhibits a 
stron g radial dependence (see iSchinnerer et al.l 
2013 and references therein). While star forma- 
tion falls on the convex side of the spiral very 
clearly between 60"<i?<80"(2<i?<3 kpc), the 
young stars and gas coincide in the outer mate- 
rial arms. Meanwhile, there is little, if any, ev- 
idence of star formation in the central 20" (in 
the nuclear bar zone; but see below) and be- 
tween 35"<_R<60" in either Ha or 24 [im emis- 
sion[3 The radial dependence in star formation 
activity in M51 causes the measured star forma- 
tion rate per unit of gas, or star formation effi - 
ciency SFE=£ SFi? /£ sas (i.e. ILerov et al.l l2008l). 
and hence gas depletion time Td ep =SFE _1 , to vary 
by as much as a factor of 6 throughout the inner 
9 kpc. This same radial var iation is evident i n 
the SFE profile calcu lated bv ILerov et~aTl (|2009| ). 
iKnapen et al.l (|1992f) identified this same pattern 
and loosely interpreted this as a dependence of the 
SFE on location relative to spiral arm resonances. 

11 Both star formation tracers show a strong deficit in the 
zone 35"<i?<60", in particular, ruling out the extinction 
of optical photons as the source of the radial variation in 
the inferred SFR. 



In Figure [TJ we compare the molecular gas de- 
pletion time Tdep to the profile of present-day grav- 
itational torques. To calculate Td ep , we divide the 
radial profile of the molecular gas surface den- 
sity T,H2 by the star formation rate surface den- 
sity Y^sfr profileo By using azimuthal averages 
of T*sfr and E# 2 we minimize the unavoidable 
bias introduced by the spatial offset between gas 
and young stars predicted by density wave theory- 
and well-known to appear along the spiral arms in 
M51. 

It is clear from Figure [T] that there is a strong 
link - shown as an inverse correlation - between 
Tdep and the gravitational torques, underlining a 
strong dependence of star formation on dynamical 
environment. Radially inflowing gas appears to be 
less efficiently forming stars. More generally, gas 
sitting near CR forms stars more efficiently than 
gas that is away from CR and in motion relative to 
the background potential. We expect this to ap- 
ply in other systems where radial outflow motions 
could potentially coincide with low star formation 
efficiency. 

In contrast to expect ations for a univ ersal 
gas depletion tim e (i.e. iBigiel et~aLl 120081 and 
iLerov etldl [20 09). it therefore appears that the 
Tdep in M51 shows a genuine, and non-monotonic, 
radial dep endence. This was pr eviously demon- 
strated by [Schuster ct al. (2007J) out to 12 kpc, 
although with a slightly different radial trend and 
using a different star formation tracer. In particu- 
lar, they find a depression in the ratio TiH2/^sfr 
in the zone of the nuclear bar. However, several 
factors complicate reliable estimation of Y>sfr, 
and thus Td ep , in this area, in particular. The 
diffuse FIR contribution from dust heating by an 
older population of bulge stars, for one, may lead 
to overestimation of the SFR. Age and/or metal- 
licity gradients in the underlying old stellar popu- 
lation also lead to uncertainties in the stellar con- 
tinuum subtraction of narrow-band Ha imaging. 

In this work, we do not find the same rise i n 
^sfr toward the center as lSchuster et al.l (|2007| ). 
possibly because any diffuse component present 
in the 24 \xm image we use is surreptitiously 
minimized in combination with the continuum- 
subtracted Ha, which shows a central deficit 
(jSchinnerer et al.ll2013[ ). In light of these uncer- 
tainties in YjsfRi we consider our estimates for 
^SFR and Td e p to be least reliable in the center. 
We note, though, that an increase in Tdep toward 
the center, such as we find, might be expected 
given the high shear characteristic of these radii 
(see Figured]). 

Throughout the remainder of the paper we focus 
primarily on the the spiral arms, which show per- 
sistent variation in Tdep that is more arguably real. 
The gap in star formation between 35"<i?<60" is 

12 The gas depletion time denned in this way is the time 
to deplete the molecular material, rather than the total 
gas, although the two timescales should be similar; the 
atomic gas surface density is a factor of 10 lower than the 
molecular gas surface density across the PAWS field of view 
(Schuster ct al. 2007;; ILerov et aLHSOOH) 
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Fig. 3. — Molecular gas surface densities throughout the 
PAWS field. Catalogued GMCs (sec Colo mbo et alJ2013al) 
are shown as gray dots, together with t he average sur- 
facc density in the arm region defined as in Colombo et al. 
(2013a) (black thick curve) and the average inter-arm 
surface density, calculated via inverse-intensity weighting 
(dashed thin curve). 

evident in multiple tracers, in contrast with the 
high Yj gas in this zone. This contrast is particu- 
larly obvious at the high resol ution of the PAWS 
observations (jSchinnerer et a,l.ll2013f ). But even at 
lower resolution, it is quantitatively clear that the 
star formation rate is lower than expected for the 
relatively high gas column; fewer stars are formed 
than expected for a global star formation law, or 
universal Tdep- Furthermore, the pattern in Tdep 
can not be attributed to variation in Xco- There 
is little to no indication that Xco varies across 
the PAWS field of view (indeed, M51 has a negli- 
gible metallicity gradient ), as re cently discussed in 
detail by Colombo et al. 2013a). Instead, we pro- 
pose that the radial dependence in Tdep reflects 
the influence of dynamics on cloud properties (see 
Section 14. 2p . following the large-scale organiza- 
tion of the ISM by the non-axisymmetric structure 
present in the disk. 

4. CMC EVOLUTION IN M51: THE 
IMPORTANCE OF DYNAMICAL 
PRESSURE ON CLOUD STABILITY AND 
STAR FORMATION 

4.1. Dynamical Pressure: Qualitative Remarks 

In the previous Section, we showed that there 
is a strong link between zones of strong gas flow 
(which, in M51, correspond to regions of radial 
inflow and negative torque) and regions where 
molecular gas appears to be inefficient at forming 
stars. This is especially striking in the inner spiral 
arms spanning 35"<i?<60", where the molecular 
gas reaches relatively high surface densities but is 
effectively devoid of high-mass star formation ac- 
tivity, in stark contrast to models that propose 
that the surface density of SFR and gas are well- 
correlated. 

In particular, Figure [2] shows that the molecular 
gas exhibits a complex response to M51's under- 



lying non-axisymmetric potential. The rates at 
which gas flows along and across the arms depend 
on position along the arm, and these rates vary 
in magnitude relative to one another depending 
on distance from corotation. The length of time 
spent in the potential well of the arm depends on 
the perpendicular component of the velocity v±, 
which is well-approximated by the radial stream- 
ing velocity v r given the tightness of the spiral. 
How far gas travels along the arm before exiting 
on the downstream side depends on the azimuthal 
streaming velocity v$. The smaller this rate, the 
more quickly gas moves in to and out of the spi- 
ral, leaving little time for clouds to acquire mass 
(i.e. through collisions/agglomeration in the spi- 
ral arm), become gravitational unstable and form 
stars. At corotation, both velocities are zero and 
gas is stationary in the rotating frame. 

We can begin to appreciate the effect of such 
motion on clouds by considering the ratio v 2 s /u 2 , 
which is a ratio of timescales: the velocity disper- 
sion a sets the 'equalization' time, or the timescale 
for clouds to reach pressure equilibrium with their 
surroundings, and the streaming velocity vs sets 
a travel, or relocation, time. In this case, Vg/Aa 2 
is a measure of the square of the ratio between the 
time it takes to pressure-equalize and the time to 
cross two cloud lengths^ With the introduction 
of non-circular motions v>a, clouds equalize less 
quickly simply because they do not remain in any 
one environment for long enough. 

This fact is underlined in the bottom panel of 
Figure [131 which shows that the ratio of equaliza- 
tion to relocation timescales is large everywhere 
except near CR where, by definition, gas is mo- 
tionless in the rotating frame. We therefore in- 
terpret the close (if qualitative) relation between 
Tdep and Vg/a 2 to suggest that moving clouds 
have insufficient time to collapse and form many 
stars. This can be viewed as an effective stabiliza- 
tion, quantifiable in terms of dynamical pressure 
Pdyn = pVg. As described in the next section, dy- 
namical pressure reduces the internal pressure of 
the streaming medium in which clouds are embed- 
ded. This reduces the surface pressure on clouds, 
countering the gravity of the cloud. 

According to the timescale arguments above, 
the dynamical pressure does not 'dissipate' in 
a cloud-crossing time, implying that dynamical 
pressure can not be ignored. Of course, whether or 
not clouds are sensitive to changes in surface pres- 
sure depends on whether clouds are dynamically 
decoupled from their surroundings. This depends 
on the internal pressure of the cloud in relation to 
the external pressure. Observations of clouds at 
high density contrast to the surrounding medium 
are usually taken to imply an extreme pressure 
imbalance, and thus little sensitivity to environ- 
ment (e.g. Krumholz. Dekel fc McKee 112012] and 
references therein) . However, Hughes et al. (2012) 
have recently shown that the internal and exter- 

13 We adopt the two-cloud length crossing time for this 
particular ratio for consistency with the formalism in £14.21 



nal pressures of clouds in M33, the LMC and M51 
(among others) track each other, even if they may 
no t be equal, similar to the cl ouds in M64 studied 
bv lRosolowskv fc Blitzl (pOOl . 

We thus have reason to expect that pressure 
plays an especially important role in M51, not 
least because the high molecular to atomic gas ra- 
tio si gnifies a very hig h midplane pressure (e.g. 
Blitz & Rosolowsky 2004). This was discussed by 
Hitschfcld ct al. (2009) who find that the hydro- 
static pressure is tightly correlated with the frac- 
tion of molecular to atomic material, indicating 
that this physical parameter determines the for- 
mation of molecular gas from the atomic phase. 

As shown in Figure [3] the surface density in 
the ambient ISM is comparable to, and in some 
places exceeds, the surface density of individual 
GMCs. This implies that the internal and ex- 
ternal cloud pressures are similar, and so we ex- 
pect that clouds embedded in this medium will 
be to some degree pressure supp orted (see, e.g., 
iKrumholz. McKee fc Tumlinson M2009f ). A large 
fraction of clouds from across the PAWS field do, 
in fact, appear to be supported by pressure ac- 
cording to their measu red virial parameters a>l 
([Colombo et al.l [2013a). Clouds from the inner 3 
kpc where the surface densities in the arm lie well 
above S// 2 =100 Mq pc -2 preferentially exhibit 
even higher a>2 than clouds from elsewhere in 
the disk. These clouds, which are comparable to 
clouds under starburst conditions, will therefore 
be keenly aware of changes in surface pressure. 

Even in lower surface density environments, 
changes in surface pressure will effect cloud equi- 
librium, making the difference between stability 
and collapse. This may be the case in other spiral 
galaxies with lower molecular gas fractions than 
M51 since, even in these galaxies, arm surface 
densities will often lie near T,qmc, even if the av- 
erage surface density does not, preventing clouds 
from completely decoupling from their surround- 
ings. For example, if we write the surface density 
of a gas disk hosting a two-armed spiral as 



S(i?,0) = S o (i?)cos(0) 2p 



(1) 



where higher (even) 2p describes a thinner spiral, 
then the average gas surface density 

1 r°° (2-nV 

< E(i?) >= — / S(i?, 0W = MR)j^l 
2tt J APp\p\ 

(2) 
will be ^0.25-0.5 times the peak surface density in 
the arm. According to the surface densities pro- 
filed by Leroy et al. (2008) we expect that, within 
the molecule-dominated zones in spiral galaxies, 
the gaseous spiral arms will supply non-negligible 
pressure support for GMCs. 

4.2. Dynamical Pressure: A Model Based on the 
Bernoulli Principle 

In the previous section we established a quali- 
tative correlation between the gas depletion time 
and dynamical pressure, which suggests that 



changes in the cloud surface pressure induced by 
gas kinematics influence cloud stability, and there- 
fore the ability of a cloud to form stars. To de- 
scribe how these changes emerge we next develop 
a simple model based on the Bernoulli principle, 
which equates an increase in gas velocity with re- 
duced pressure. 

Specifically, we work under the assumption of 
isentropic flow of a compressible isothermal ideal 
gas and relate gas properties along the stream- 
line between two zones using the compressible 
Bernoulli equation. This is a simplification and 
ignores the fact that gas in reality undergoes 
shocks and other forms of energy loss; gas in per- 
turbed disks moving at non-zero velocity relative 
to the background rotating frame will exhibit tran- 
sonic speeds (Mach number Jl/>1). The shock 
can be weakened by strong gas self-grav ity (e.g. 
IKim. Kim fc Ostrikerl [2001 IWadal 12001 and we 
expect this to be true in M51. But in the frame of 
the shock, entropy is still not conserved and the 
values for the quantities assumed in the expres- 
sions in the following section will depart from the 
values assumed (the less so the weaker the shock). 

In particular, under the assumption that the 
dominant source of pressure is turbulent in nature 
so that P=P&turbi we assume that the internal tur- 
bulent motion of the gas is uniform throughout 
the disk and take the average vel ocity dispersion 
<7«nr-h~ 10 kms -1 observed in M51 (jColombo et al.1 
l2013bf) as representative of this motion. This value 
takes into account the contribution of unresolved 
bulk motions to the observed linewidth, which 
we subtract (f ollowing the formalism adopted by 
iPety et aLll2013[ ). Removing these motions yields 
a much more uniform velocity dispersion, although 
there is still a factor of 2 difference in the turbulent 
motions between the arm and interarm, likely be- 
cause of the shock. Since this factor of 2 contrast 
is smaller than the density contrast between arm 
and interarm (typically 5-10), we ignore it. We 
can more closely approximate reality by letting 
the density depend on pressure, i.e. by assuming 
compressibility. 

Only with a better observational handle on en- 
ergy losses and shocks in gas (which is minimal 
at present) will it be possible to invoke an expres- 
sion equivalent to Bernoulli but with the change in 
entropy included. This is something for which nu- 
merical calculation is better equipped and so we 
proceed with our simple calculation below with 
the aim of providing at least a coarse description 
of the phenomenon. 

4.2.1. The stable GMC mass threshold 

In this section we compare the maximum cloud 
mass stable against collapse in two different flows, 
i.e. in the arm and inter-arm regions. Gas mov- 
ing at high relative velocity, tunneled along and 
through the arms, will result in a reduced ex- 
ternal surface pressure for clouds forming in the 
arm. This then increases the maximum mass 
stable against collapse, as related by expression 
for the Bonnor-Ebert mass MbeO(. P -1 / 2 calcu- 



lated from pressure, rather than virial, equilib- 
rium, i.e. including a (thermal) surface pressure 
term. FiFi As a result of the dynamical pressure 
1) any collapse-unstable clouds that form will be 
of higher mass than virialized clouds and 2) fewer 
clouds will be collapse-unstable, leading to lower 
SFE. 

We can estimate the fractional increase in the 
Bonnor-Ebert mass as the result of dynamical 
pressure by equating the energy densities of gas 
in two zones with and without strong streaming 
motions (e.g. in the arm and interarm) in the non- 
intertial (rotating) reference frame. For simplicity, 
we use the compressible Bernoulli equation to re- 
late gas properties along the streamline between 
the two zones. For pressure P dominated by tur- 
bulent motion with dispersion a (i.e. Pi=Pi<J ; 
note that this ignores an enhancement in cloud 
internal turbulence during passage through the 
arm) , the pressure ratio of gas in zone 1 and zone 
2 can be expressed as 

ln p;=2^"2^ + ^2-(^2-Si/ l i) (3) 

with hi the gas scale height and Sj the gas sur- 
face density in the two zones. The first two terms 
on the right are the dynamical pressures, while 
the third represents the change in the hydrostatic 
pressure from zone to zone. Omitted terms in- 
clude the stellar density factor in the hydrostatic 
pressure and other energy losses/gains. 

As we are near enough to hydrostatic equilib- 
rium so that Ei/E2~/i2//ii, the third term on the 
right is negligible and the first two terms domi- 
nate. 

Since fractional increase in the stable (Bonnor- 
Ebert) cloud mass is 
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From this point we can see that streaming mo- 
tions raise the cloud stable mass above the virial 
mass by 



log 
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4ct 2 In 10 ' 



(5) 



where the non-circular streaming in zone 2 is given 
by vs and zone 1 is identical to zone 2 except now 
its clouds are free of streaming (and dynamical 
pressure) and thus assumed to be virialized. 

Within gas rotating in the plane of the galaxy, 
the stable mass is raised for all clouds under- 

14 Note that the Bonnor-Ebert mass is equivalent to the 
Jeans mass ignoring any external or dynamical pressure. 

15 Even if the surface term were of kinetic origin (rather 
than thermal, as assumed in the expression for the Bonnor- 
Ebert mass) we would expect an analogous change in stable 
m ass with a change in surface press ure, such as considered 
bv lBallesteros-Paredes et al.l lf20T3 ) 



TABLE 1 
Environment-dependent correlations 



environment 


7 


Tdepfi 


x 2 


I 


-1.75±0.16 


0.15±0.06 


5.30 


II 


-1.1±0.18 


0.73±0.21 


0.88 


III 


-1.54±0.24 


0.89±0.28 


1.98 


IV 


-1.07±0.05 


0.72±0.13 


6.85 


V 


-1.33±0.08 


0.14±0.04 


7.43 


environment 


7 


Tdepfi 


x 2 


I 


-1.82±0.33 


0.52±0.12 


21.81 


II 


-1.08±0.19 


0.79±0.15 


0.96 


III 


-1.48±0.22 


1.54±0.25 


2.16 


IV 


-1.07±0.05 


0.75±0.11 


6.6 


V 


-1.15±0.04 


0.43±0.06 


13.35 


environment 


7 


Tdep.O 


x 2 


I 




0.11±0.02 


9.25 


II 


. . . 


0.52±0.06 


26.74 


III 


. . . 


0.68±0.11 


14.03 


IV 




0.14±0.02 


112.75 


V 




0.04±0.01 


108.7 



NOTE. -— Estimates of the cloud mass spectrum 
index 7 and fiducial gas depletion time Tdepfi mea- 
sured from the slope and intercept, respectively, of 
the best-fitting straight line in each of the 5 distinct 
dynamical environments indicated in Figure [4] The 
X 2 value of the fit is given in the right-most column. 
The five environments are, in ascending order, nuclear 
bar, molecular ring, spiral inside CR, spiral outside 
CR, and outer material pattern. The top third of the 
table lists the results of the fits shown in Fig. [4] The 
middle third lists the results of fits to similar points, 
but that take into account inter-arm streaming, and 
the bottom third lists the Td ep fi measured assuming a 
fixed slope corresponding to 7=-1.72, the average in- 
dex measured via direct fitting of the cloud mass spec- 
tra in these same environments (Hughes et al. 2012). 



going non-circular motion. The stable mass for 
clouds moving at streaming velocity i>,g^l5-20 km 
s , typical of spiral arm streaming in M51, will 
be at least 1.5-2 times higher than for vital- 
ized clouds with similar surface densities and line 
widths cr<10 km/s. As a result of the reduced sur- 
face pressure, collapse-unstable clouds that go on 
to form stars will be on average 1.5-2 times more 
massive than virialized clouds. 

The difference from arm to intcrarm cloud 
masses should be similar, given the weaker stream- 
ing in the interarm. Conservatively allowing a fac- 
tor of 1.5 difference in vs, then we expect bound 
cloud masses to be on average 0.5 dex higher in the 
arm than in the interarm, very nearly what is mea- 
sured from th e cloud-decomposed P AWS CO (1-0) 
emission (see lColombo et aLll2013al : evident from 
the horizontal offset in the mass spectra for clouds 
in the two zones). Again, this ignores real varia- 
tion in the gas dispersion from arm to interarm, 
but we note that the change in measured cloud 
line- width betw een these two environm ents in M51 
is only modest (jColombo et alJl2013af ) 



4.2.2. Impact on Tdep 



In the previous section we showed that, as a 
consequence of strong non-circular streaming mo- 
tions, the collapse-unstable clouds that go on to 
form stars will be on average 1.5-2 times more 
massive than virialized clouds. This suggests that 
star-forming clouds in the arm (where streaming 
is high) preferentially host massive star formation. 
But at the same time, any clouds in the arm be- 
low ~2 times the virial mass will not collapse and 
form stars. This naturally translates in to a de- 
crease in the star formation efficiency, considering 
the change to the fraction of the total gas reservoir 
in the form of collapse-unstable clouds per free fall 
time. 

Specifically, assuming that the process (or pro- 
cesses) responsible for the build-up of GMCs re- 
sults in a mass spectrum 



N(> M) oc 



/ M 



\M„ 



7+1 



(6) 



then from the mass in (unstable) clouds above the 
stable threshold 

M(> M stable ) = f """ ^MlhldM (7) 



M 3tM e 



dM 



we see that the mass fraction of such clouds F = 
M(> M stable) /Mtotal in a region of total mass 
M total will change as 



dF=(l + 7) 



/ M 



\M n 



dM 

Mtotal 



(8) 



with a change to the stable mass by an amount 
dM. 

Now defining the average unstable cloud mass 
<M C >=M{> M stab ie)/N, then 



F v 



■7) 



dM 



<M C > 



(9) 



We can cast this in terms of the SFE by rec- 
ognizing that the collapse unstable clouds-those 
above the stable threshold-are the clouds respon- 
sible for star formation. Specifically, we write 
SFE=F/3 where ji is the fraction of a cloud that 
is converted in to stars per free fall time (typ- 
ically less than 5% ; e.g. iKrumholz fcMcKeel 
120051: lMurravll20TlT ). In this case, the fractional 
change to Td ep , the inverse of SFE, is given by 
dTdep/rdep = —dF/F whereby 



In 



Tdep 
Tdepfi 



-(1 + 7) In 



M 



<M r > 



-(1+7) In 



= -(l+T>S 



M 

(10) 



using that <M c >~M V i r in the absence of stream- 
ing and where constant of integration Tdepfi is the 
unperturbed gas depletion time. By design, we 
take Tdepfi to be the depletion time in the case of 
virialized clouds. 
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Fig. 4. — (Top Left) Scatter plot of In Tdep vs. Vg/ia extracted in 2.4" bins from the curves at right. The symbols and 
color scale varies as at right, according to dynamical environment. Best-fit lines for each environment are overlayed. (Top 
Right) Profiles of ln-r^p (solid black) and i>|/4<r 2 (dashed black) throughout M51. The measurements of vg and a used 
here are shown in the bott om panel of Fig u re [2] E rrors have been omitted for clarity. Dynamical environments defined 
in § IB. 21 (and described by Colombo ct al. (2013a)) are designated by color overlays: nuclear bar (red), molecular ring 
(orange), spiral inside CR (green), spiral outside CR (blue), material pattern (purple). (Bottom) Scatter plots of In Td ep 
vs. Vg/Ao 2 separated by environment. 



4.2.3. Comparing Model Predictions to 
Observations 

According to eq. (fTUj) . the gas depletion time, or 
inverse SFE, along the arm should proportionately 
follow the raised stable mass threshold in the pres- 
ence of spiral streaming as measured by Vg/Aa 2 . 
These two quantities arc plotted as a function of 
radius in the left panel of Figure IU using our mea- 
surements of the magnitude of streaming motions 
shown in the bottom panel of Figure [2] and the 
observed velocity dispersion cr|3 which is found 
to be relatively uniform across the disk, near 10 
km s _1 (see Figure [5]). With our very simple pre- 
diction we are able to reproduce the observed pat- 
tern of star formation; the radial variation in the 
Tdep nearly echoes the radial dependence in the 
kinematic term v 2 s /a 2 to within a radially varying 
scale factor -(7+I). 

The close match (albeit with scatter) is under- 
lined on the right panel of Figure 21 where points 
fall along environment-unique lines with slopes - 

16 That is, we assume <r describes the overall state of the 
gas (i.e. the combined thermal and kinetic pressure on the 
gas), including the contribution from turbulent motions. 
Note that, independent of streaming, gas with higher <r is 
expected to contain clouds with higher virial masses. 



(I+7), according to eq. [KjJ Table Q] lists the 
slopes and intercepts of the best-fitting linear rela- 
tionship in each of the five environments indicated 
in Figure 0] determined with the FITEXY routine 
(|Press et al. 1992) assuming uncorrelated errors in 
both variables. We adopt 0. 2 dex uncert a inty i n 
the measured Tdep, following iLerov et al.l J2008), 
and the errors plotted in Figure [2] for the uncer- 
tainty in v^/Aa 2 . Despite the presence of scatter, 
the fitted relations are robust to changes in the ra- 
dial range over which each dynamical environment 
extends. Increasing or decreasing the size of each 
by 10% leads to typically less than 10% variation 
in the fitted slope and intercept. 

The average index 7=-1.32±0.26 measured from 
the best-fit slopes is consistent with (if slightly 
lower than) the i ndices measured in th e same set of 
environments bv lHughes et al.ll2013al via direct fit- 
ting of the cloud mass spectra (<7d/>=-1.72±0.39 
on average). A 2° error in the adopted inclina- 
tion angle can easily explain this difference (see § 

MM- 

Moreover, the measured indices are remarkably 
compatible with the genuine, if modest, trend in 
the 'direct-fit' jdf with dynamical en vironment. 
In the study of iColombo et all (|2013af) variations 
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Fig. 5. — (Top) Cloud mass spectrum power-law index 
from within the series of M51 dynamical environments des- 
ignated in Figure [4] The center of each environment is 
taken as the radial position of each measurement. Indices 
estimated from two independent approaches are shown: 1) 
the range spanned by the minimum and maximum indices 
measured from direct fits to the cloud mass spectra in each 
environment is show n in light gray (Colombo ct al. 2013a; 
Huehcs ct al. 2013a). The gray horizontal line marks the 
average 7=-1.72. 2) The indices implied by the best-fit 
slopes in Figure [4] in each environment are shown in solid 
black with errors given by the statistical error in each 
measured slope. Best-fit slopes accounting for inter-arm 
streaming velocities are shown in thinner dashed line. (Bot- 
tom) Fiducial gas depletion time Td ep q in the same series 
of M51 dynamical environments measured from the best- 
fit intercepts in the right panel of Figure [4] (solid black), 
the best-fit intercepts when interarm streaming velocities 
are included (thin dashed black) and the best-fitting Td ePi o 
with fixed 7=- 1.72 (dotted gray). 



in mass spectrum shape between the center, arm 
and interarm environments is interpreted in terms 
of cloud formation, growth and destruction mech- 
anisms. We retrieve the record of these processes 
across the environments sampled here with our 
simple expression (but see § I4.2.4J) . 

Meanwhile, the best-fit intercepts of each 
environment-dependent relation correspond to a 
fiducial gas depletion time Td ePl o = 0-5 Gyr, on av- 
erage. This value represents the time to con- 
sume the molecular gas in the absence of spiral 
streaming and is notably shorter than the aver- 
age depletion time across the PAWS field of view, 
<Td ep >= 1.5 Gyr. The latter value is consistent 
with molecular gas depletion times measured in 
the THINGS sample of nearby galaxies by Bigiel 
et al. (2011) and Leroy et al. (2008, 2012). In 
some environments the fiducial Td ep ,o appears to 
be even shorter than 0.5 Gyr (see Figure [5] but 
see ? I4.2.4[) . Later in §[5] we develop a picture that 



attributes departures from a fiducial, universal gas 
depletion time to motions in non-axisymmetric 
potentials. 

We emphasize that the trend traced out in the 
right panel of Figure H] does not arise with ra- 
dial gradients in the kinematic term v^/Aa 2 or 
Tdep across the PAWS field. Variation appears 
on smaller (arcsecond) scales, and the agreement 
between predicted and observed Tdep is genuine. 
This may be surprising, given that our predic- 
tion is oversimplified, not least because it implic- 
itly assumes a fixed, glo bal power-law m a ss spec - 
trum. As revealed by iColombo et al.l (|2013aD . 
cloud mass spectra in M51 show marked devi- 
ations from power-law behavior and uniformity 
along the arms coincident with changes in the 
dynamical environment defined here. These dif- 
ferences appear to indicate real variations in the 
mechanisms behind GMC formation, but they 
may also reveal how dynamical pressure plays a 
role in shaping/modifying an underlying global 
power-law mass spectrum (and lognormal inten- 
sity pdf; see Hughes et al. 2012). Detailed mod- 
elling will be required to better understand and 
distinguish these scenarios. 

4.2.4. Sources of Uncertainty: Sensitivity to 
Observables 

Several factors influence the tightness of the cor- 
relation between Tdep and Vg. In the previous sec- 
tion we demonstrated that the overall scatter is 
dominated by real environmental variation in 7 
and Tde P ,o- Within a given environment, the scat- 
ter that emerges tends to be on the order of 5- 
7%, according to the errors on the measured slope 
and intercept of the best-fit linear relation. Both 
within and among environments, the correlation 
is subject to observational uncertainty, namely to 
the accuracy in the star formation and gas den- 
sity tracers, the kinematic parameters assumed 
to deproject the line-of-sight velocities and the 
pitch angle adopted in the calculation of the non- 
circular streaming motions. 

The reported statistical errors in the fitted Td ep ,o 
and 7 primarily reflect uncertainty in the assumed 
pitch angle and rotation curve, which together de- 
fine the errors in v 2 s used as weights in the reduced- 
X 2 fit to the linear relation in eq. (JTDJ). In testing 
we found that pitch angle variation is responsible 
for the largest uncertainty in the decomposition 
of the line-of-sight velocity field: by comparison, 
the radial and azimuthal streaming components 
are surprisingly robust to changes in the kinematic 
parameters, i.e. to the major axis position angle 
and inclination. We estimate that PA variation of 
±5° introduces 10% uncertainty in the magnitude 
of the streaming motions and therefore up to 20% 
change in the best-fit slope. Velocities extracted 
via solution in the spiral arm frame exhibit the 
expected sini dependence. We expect inclination 
uncertainty, unlike pitch angle, to apply globally 
(under the assumption of an unwarped disk), re- 
sulting in an overall shift in the fitted Td e p,o and 7. 
A slightly higher inclination of i=23° would bring 
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the average 7 in to perfect agreement the average 
index measured by Hughes et al. (2012). 

To construct our estimate of the star formation 
rate, we use a combination of Ha and 24 fim emis- 
sion. This accounts for obscured and unobscured 
star formation but does not explicitly include a 
correction for the diffuse component of the 24 /im 
emission (e.g. as explored by Leroy et al. 2011), 
which is thought to arise with dust heating by an 
underlying older population of stars. The central 
bulge region (but also the spiral arms) most likely 
incorporates some level of dust emission that does 
not trace young stars, in light of the enhanced old 
stellar surface densities in this zone (i.e. as demon- 
strated in the bulge region of M31 by Groves et 
al. 2012). Underestimation of Tdep in the center 
as a result is likely responsible for the particularly 
steep power law index estimated for the nuclear 
bar environment, as well as the low Td e p,o there. 

A low Td e p m the center could alternatively arise 
from underestimation of the molecular gas surface 
density. Given the excellent sensitivity to low- 
level CO emission in the PAW S data set and little 
expected variation in Xco (|Hughes et al.ll2013al 
iColombo et a l. 2013a) the latter would most likely 
occur only if the CO fails to trace a warmer molec- 
ular phase, e.g., as more regularly observed in the 
centers of galaxies (Dumas et al.). The gas de- 
pletion time in the outermost material arm envi- 
ronment may also be artificially lowered due to 
the increased contribution from atomic gas com- 
ponent omitted at these radii (although this is still 
below 10% at the edge of the PAWS field; see the 
profiles in Leroy et al. 2008). 

Finally, the quality of the correlation depends 
on the reliability of our assumptions. We ig- 
nore any non-uniformity in the turbulent motions 
from arm to interarm as well as energy losses 
that occur in the spiral shock. Our isothermal 
assumption further ignores complex heating and 
cooling within the ISM. In reality, clouds will 
show variation in temperature between, e.g., the 
spiral arm and interarm regions. New simula- 
tions that implement heating and cooling pre- 
scriptions find dispersions in the range 3-7 km/ s 
([Hopkins. Quataert fc Murravl2012l : lDobbsl201lT ). 
According to eq. (fT0|). an error of this magnitude 
(~50%) would correspond to an uncertainty of 0.1 
dex in the gas depletion time (25% uncertainty in 
Tdep)- As this i s the typica l mea surement uncer- 
tainty in Tdep (|Lerov et al.l 120081 ) our formalism 
can conversely accommodate 50% uncertainty in 
the velocity dispersion. Given the current data 
quality, we feel that the isothermal assumption 
is valid for our purpose, although future obser- 
vations (i.e. with ALMA) will make it possible to 
test such a scenario. 

Numerical calculations are clearly necessary to 
perform a detailed energy balance and follow the 
evolution and equilibrium of individual clouds 
during passage through the spiral arm. Still, the 
fact that we observe a correlation at all provides 
some indication that our simple expression is good 
to first order, even neglecting shocks and energy 



losses. 

The form of the correlation we consider here 
(plotting In Tdep vs. u|/4cr 2 ) by design represents 
the impact of dynamical pressure on the length- 
ening of the gas depletion time with respect to 
that of virialized clouds. But this ignores that, 
upon entering the arm, some clouds may them- 
selves not be in virial equilibrium, i.e. due to non- 
zero streaming in the interarm. Even accounting 
for non-zero interarm streaming introduces little 
change to the correlation from environment to en- 
vironment. We recalculated all fits including the 
inter-arm streaming velocities, i.e. according to 
eq. [4) where now the lengthening in gas deple- 
tion time in the arm is relative to the depletion 
time in the interarm. The 7 and Tdep,o associ- 
ated with the best-fit slopes and intercepts in each 
of the five environments under consideration are 
listed in the middle part of Table [TJ We find on 
average Td e p,o=0-66±0.26 and 7=-1.56±0.74 and 
a trend with environment that similarly resem- 
bles that obtained from direct fits to the mass 
spectra, although with slight exaggeration (see the 
blue points in the top and bottom panels in Fig- 
ure [5]). Given that the interarm velocities are 
smaller than the arm velocities, this adjustment 
is responsible for less deviation from the simplest 
case than when we impose a global power law mass 
spectrum. Estimates for Td ep ,o accounting for the 
inter-arm streaming are more similar to the val- 
ues estimated including only arm streaming than 
to those where the slope is fixed to 7=-1.72 (see 
the entries in the bottom part of Table [T] and plot- 
ted in Figure [5]). 

5. DISCUSSION 

The close link between the predicted stable 
cloud mass and Tdep described here strongly im- 
plicates the influence of dynamical pressure on the 
formation and evolution of the GMCs observed in 
M51. Indirect evidence for the role of surface pres- 
sure in cloud equilibrium is also seen in the lack 
of a clear size-line width relation for cloud s in the 
arm (and interarm; IColombo et al.ll2013at Hughes 
et al. 2012b), which, when observed, is taken as 
evidence for viriali zed clouds. Our fin dings reiter- 
ate the warnings of lShettve t al. (2011) against the 
use of the virial parameter for determining cloud 
properties and quantifying boundedness. The ex- 
ternal pressure-and the decrease in surface pres- 
sure in the presence of strong streaming motions 
in particular-must be taken into account. 

The cloud-decomposed emission in the PAWS 
data cube shows other consistencies with the pre- 
dictions for the impa ct of dynamical pressure. 
IColombo et al.l ()2013aD present evidence that the 
GMCs in the arms are more massive than those 
in the interarm by about an order of magnitude, 
clear by the offset in the mass spectra for clouds in 
the two zones. Meanwhile, a larger fraction of the 
inter-arm g as is in the form of cl ouds than in the 
spiral arm (|Colombo et al.ll2~013a|) . In our picture, 
regions with high streaming velocities should have 
fewer collapsing clouds, and the few that exist will 
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be of higher mass than clouds that are virialized. 
At present, the tight link between gas flows 
and star formation is based on only one galaxy, 
and tests of similar systems will be necessary to 
firmly establish the role of dynamical pressure. 
In this discussion section, we compare dynamical 
pressure with other sources of cloud stability and 
highlight several aspects of star formation that 
can be explained uniquely by our picture. 



5.1. Other potential sources of cloud stability 

In the previous section we revealed an 
environment-dependent correlation between gas 
depletion time and dynamical pressure. This 
arises from the dependence of the star formation 
efficiency on the mass fraction of collapse unsta- 
ble clouds per free fall time, which varies with lo- 
cation in the disk; the mass fraction of collapse 
unstable clouds is specified by the environment- 
varying cloud mass spectrum index, together with 
the stable mass threshold set by the dynamical 
pressure. In our picture, dynamical pressure ef- 
fectively stabilizes the clouds and thus prevents 
all but those above the raised stable mass thresh- 
old from collapsing to form stars. But there are 
several other potential sources of cloud stability. 
We review these here and consider the possible 
role shocks might play in stimulating the observed 
pattern of star formation. We conclude that none 
of these cases can adequately explain the observed 
pattern in the SFE or gas depletion time, leaving 
dynamical pressure as the most compelling source 
of stability. 

5.1.1. Shear 

As the locations of strong non-circular mo- 
tions and enhanced surface density, spiral arms 
are prefe rentially sites of reduced shear (as ar- 
gue d by Elmcgrccn 1987; lElmegreenl 119941 an d 
see IKim fc Ostrikerlbooa lkim fe Ostrikerl 120061) . 
This also makes them favored sites of star forma- 
tion since, in the absence of shear to counter self- 
gravity, gas can form coherent cloud structures 
that become gravitationally unstable. 

The link between enhanced shea r and r educed 
star for mation has been stud ied by ISeigad (12 005) 
(and see iWeidner et al.|[20lol and iHocuk fc Spaansl 
1201 If) . While this may hold on global scales, 
Figure [6] suggests otherwise for M51. The link 
between reduced shear and star formation does 
not hold uniformly along the spiral arms on two 
counts: 1) Shear as measured by the Oort param- 
eter A =l/2(v l f ) /R - dv^/dR) in the spiral arm is 
very clearly non-zero along the length of the arm. 
Shear in the spiral arm only counters the local 
background shear (from the rising rotation curve; 
see O in Figure [12]) between 18"<i?<28"and fur- 
ther out near i?~85"and R^lOO". In these outer 
regions, in particular, it appears that shear is re- 
sponsible for the lengthening of the depletion time. 
2) The region 35"<i?<60" where we find a pro- 
nounced increase in Tdep (or decrease in SFE) is 
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Fig. 6. — (Top) Profiles of the background shear due to 
differential rotation (gray dashed line) and the shear in the 
bar and spiral arm regions, calculated from our measure 
of v^ in Figure [2] (black curve) . Thin black lines represent 
the rms dispersion in solutions with ±5° variation in the 
assumed pitch angle i p . (Bottom) For reference, the radial 
trend in the molecular gas depletion time is depicted in 
solid black, repeated from Figure [T] 



characterized by a similar, if not lower, degree of 
shear as in the neighboring zone at larger galacto- 
centric radius, which is characterized by relatively 
more star formation. Moreover, as examined in 
Appendix [D] the threshold for cloud instability 
and collapse in the presence of shear (as well as 
Coriolis and tidal forces) is exceeded by the ob- 
served gas surface density everywhere along the 
length of the arms. To prevent star formation on 
cloud scales between 35"<i?<60" another source 

of stability is clearly required; 

Our finding echos that of iDib et al.l (|2012l ) who 
considered the role of shear on the scale of individ- 
ual molecular clouds in the Galactic Ring Survey. 
In all cases, they find that shear does not balance 
gas self-gravity and conclude that turbulent mo- 
tions, stellar feedbaclfH and/or magnetic fields are 
the principle agents of cloud stability. 

5.1.2. Turbulence and Magnetic Fields 

To produce the observed pattern in Td ep in 
M51 we might similarly expect turbulent mo- 
tions to be enhanced within clouds in the zone 
35"<i?<60" compared to clouds from neighboring 
zones. Our high resolution PAWS observations, 

17 Star formation driven winds from successive star for- 
mation events are thought to deposit the energy required 
to sustain turbulence (which oth erwise dissipates quickly; 
i.e. Murrav. Quatacrt & Thompson 2010) 
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where the line-of-sight velocity profiles sample tur- 
bulent velocities above ~2 km s _1 ()Petv et al.1 
2013? iColombo et al.|[2~013al ) on cloud scales, do 
not show such a pattern, however. The veloc- 
ity dispersion shows no pronounced change along 
the arms (see Figure [5]). The virial parameter 
a, which measures the balance between the in- 
ternal kinetic energy and gravity of clouds, shows 
no significant difference between this region of the 
spiral ar m and the nei ghboring segment, beyond 
R=60" (|Colombo et al.ll2013al) . Either turbulent 
motions are not responsible for the radial variation 
in Tdep shown in Figure [1] or they are present but 
fall below our velocity resolution. Below we will 
consider shocks and their influence on cloud veloc- 
ity dispersions in more detail. Here wc note that, 
given the marked absence of star formation in this 
zone, cloud support can not come from turbulence 
driven by stellar winds/feedback. 

Magnetic fields are another potential influence 
on the organization of the ISM, but their role in 
cloud stabilization is not clear (or observationally 
established) at this point. Since they should per- 
vade the ISM and clouds equally, it seems unlikely 
that they are responsible at the cloud level for the 
trends we see here. 

5.1.3. Shocks 

Spiral arm shocks, on the other hand, have been 
long appreciated for their role in the formation 
of G MCs and the triggering of star formation 



e.g. iRoberts fc Stewartill987t iGittins et al.H2 003; 



Gittins fc Clarkd 120041 [Kim fe Ostrikerl (I2002D : 
Dobbsl 120081 ) . Shocks bring the gas surface den- 



sities high enough to initiate gravitational insta- 
bility and at the same time favor massive GMC 
formation via collision and agglomeration. Clouds 
formed as a result of shocks are predicted to hav e 
elevated velocity dispersions (|Bonnell et al. 2006), 
with clear implications for cloud boundedness and 
hence the global patterns of star formation. 

In this shock picture, variations in cloud prop- 
erties would also presumably emerge along the 
spiral arms depending on the properties of the 
shock, or as a result of variations in the balance 
between shear, self- gravity and agglomeration (as 
well as the properties of the pre- shock medium, i.e. 
iDobbs fc BonnelH 120071 : [Dobbs 2008). For exam- 
ple, the zone we identified with very little star for- 
mation might preferentially host unbound clouds 
formed via agglomeration in the shock rather than 
gravitational instability. However, there are sev- 
eral observations that counter this interpretation, 
as summarized below. 

First, the surface density in this region is among 
the highest along the arms. In comparison to the 
other zones with star formation-and which pre- 
sumably host bound clouds-this zone should like- 
wise favor cloud formation via gravitational insta- 
bility. 

At these high surface densities, unbound clouds 
might still emerge if agglomeration builds on the 
masses of clouds already formed by the dominant 
process of instability in self-gravitating gas (as re- 



alized in the simulations of IDobbsl I2008D . The 
properties of these clouds are predicted to depend 
on the strength of the shock, with stronger shocks 
leading to the a ccumulation of more massive, 
unbo und clouds (jDobbs fc Bonnelll 120071 : IDobbsl 
[2008h . 

This case could potentially mimic a varying 
gas depletion time since, from this point, stars 
formed from individual bound clouds embedded 
in a larger unbound structure do so at an over- 
all low efficiency. We might therefore expect the 
inverse of Td ep (the SFE) to be well-matched to 
the ratio of total mass in clouds to the total mass 
in the arm (i.e. a low mass of potentially star- 
forming material where Tdep is long). But Figure 
[7] shows no such behavior. Instead, the mass in 
clouds relative to the total mass seems even el- 
evated where the star formation is preferentially 
low. This is very strong evidence that the agglom- 
eration of unbound cloud associations via shocks 
is not responsible for the observed pattern in Td ep , 
especially since it conservatively assumes that all 
catalogued clouds will collapse to fo r m star s. 

As analyzed by IColombo et al.l (|2013al l. the 
shape of the cloud mass spectrum provides a 
record of the processes responsible for the forma- 
tion (and destruction) of clouds in M51. GMCs 
in the spiral arms are predominantly formed by 
gravitational instability, although agglomeration 
plays a role in organizing the ISM into cloud as- 
sociations along these structures. Notably, the 
cloud mass spectrum shows no si gnificant change 
in ove rall shape along the arm ([Colombo et al.l 
l2013al ). This suggests that the mix of instabil- 
ity and agglomeration in cloud formation is equal 
in arm environments, namely between the spiral 
segment neighbored by star formation on its con- 
vex side and the segment 30"<i?<60"with little to 
no star formation. The pattern in star formation 
therefore does not appear to be tied with cloud 
formation mechanism. 

5.1.4. Triggered star formation 

Star formation triggering by spiral-arm shock 
dissipation is also very clearly not responsible for 
the observed pattern in the SFE. Since stronger 
shocks are accompanied by larger velocity gradi- 
ents and thus dissipate more energy, we might 
have expected relatively short gas depletion times 
in exactly the locations of strongest spiral shock 
(i.e. anywhere away from corotation), unlike what 
is observed. This picture would also require the 
shock strength to vary from one segment of the 
arm to the other (i.e. 40-60" compared with 60- 
80") to explain the difference in the star forma- 
tion properties between these tw o zones, but this 
does not appear to be the case (|Schinnerer et al.l 
2013). Taking the spiral shock strength as the 
size of the offset between the molecular spiral arm 
from the potential minimum traced by the peak 
in t he old stellar light distri bution imaged at 3.6 
fj,m, iSchinnerer et al.l (|2013l ) find only a modest 
difference between the zones 40-60" and 60-80". 
We therefore conclude that, although spiral- arm 
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Fig. 7. — Histogram of the ratio of the molecular mass in 
catalogued clou ds M c i ou d s (est imated from the cloud CO 
luminosity; see Colombo ct al. 2013a) to the total (lumi- 
nous) mass in the arm M arm measured in 4.5" radial bins. 
Thin lines delineate the dispersion in cloud masses within 
each bin. 



shocks may be active in M51, it is not clear that 
they have a strong, direct influence on the trigger- 
ing of star formation. 

5.2. Implications for patterns of star formation 
in observed and simulated galactic disks 

Up to this point, simulations and semi- 
analytical models of GMC formation and evolu- 
tion rarely consider non-axisymmetric potentials, 
and so streaming motions and dynamical pressure 
have been neglected. Variations in the star for- 
mation efficiency and cloud stability have mostly 
focussed on the role of tur bulent support and 
turbu lent energy dissipation rtKrumholz fcMcKed 
2005 ; iMurray. Quataert fc Thompson! 120101: 
Fcldman n et al.l 120111 ). The few simulations with 



fixed spiral potentials have thus far focused only 
on very tightly-wound spirals with pitch angles 
«p<10° and so have not probed the regime of 
strong streaming motions. In addition, these 
simulations consider the Bernouli principle more 
for its role in spur formation via the Kelvin- 
Hclmholz instability, which has subsequently 
been found ineffective in real ga laxies compared 
to th e magneto-Jeans instability (K im fe Qstriken 
2002). 

From the observational perspective, we can al- 
ready infer the role of dynamical pressure in 
galaxy disks. The strong dynamical dependence of 
the conversion of gas in to stars presumably con- 
tributes to the frequent lack of a clear signature 
in the angular cros s correlation between gas and 
youn g star tracers (jFovle et al.ll2011t lEgusa et al.l 
2009) . The angular offset between the young stars 
and gas predicted as a result of the propagation 
of the spiral pattern (combined with the timescale 
for star formation; e.g. Roberts & Stewart 1987), 
will not be observed, even where the gas column 
is high, if the gas is dynamically prevented from 
forming stars. 

The influence of dynamical pressure is also man- 
ifest at the global level. We suspect that it keeps 



stronger grand-design spirals from forming an ex- 
cess of stars relative to weaker, flocculent spirals 
with similar global molecular gas surface densities, 
as observed (see Elmegreen & Elmegreen 1986). 
Whereas stronger spirals result in locally higher 
gas surface densities, and might therefore be ex- 
pected to exhibit higher star formation rate sur- 
face densities, they also locally stabilize clouds 
through dynamical pressure. These two effects 
presumably together keep the global star forma- 
tion rate independent of spiral strength. An in- 
creased velocity dispersion in the strongest spiral 
arms will achieve the same result, but we note that 
a distinction from dynamical pressure is difficult 
to establish, as observations at resolutions above 
the size scale of GMCs will also reflect unresolved 
streaming motions. 

The strong internal gas flows associated with 
galaxy interactions suggests that dynamical pres- 
sure may be a common feature in such scenar- 
ios. M51 and its companion are an ideal test-bed 
for exploring the impact of interactions on gas or- 
ganization and subsequent star formation. Our 
consideration of the impact of M51s companion 
galaxy has been so far indirect, concerning the 
non-axisymmetric structure across the PAWS field 
that the interaction helped to stimulate and shape. 
Features that are more directly influenced by the 
companion, namely the "bridge" that extends be- 
tween the two galaxies and the opposite "tail" , are 
beyond the edge of the PAWS field of view, pre- 
cluding the study of the ISM organization and the 
properties and progression of star formation with 
the same detail as we have undertaken at smaller 
galactocentric radius. As the ISM composition is 
changed at large radii, as well, future insight in to 
the nature of the impact of the interaction will re- 
quire extended, high resolution kinematic tracers 
of the molecular and atomic ISM phases. 

5.2.1. Relevance for barred galaxies 

The scenario we find here - low star forma- 
tion efficiency in the CO-bright inner spiral arm 
segments - can be thought of as comparable 
to the dust lanes along the length of bars typ- 
i cally observed to harbor little star formation 
(|Sheth et all2005t iComeron et al.ll2009D . Star for- 
mation is thought to be inhibited in the strongest 
bars because of str o ng shocks and h igh shear 
(lAthanassoulal 119921: iSheth et al.l 120051: but see 
Zuri ta et al.l ( 20041 ) who suggest shocks help trig- 
ger star formation). In our picture, on the other 
hand, the strong radially inward streaming mo- 
tions that accompany strong bar torquing stabi- 
lize clouds, shutting off star formation. We find 
this to be a compelling interpretation since the 
two conventional kinematic sources of cloud stabi- 
lization - enhanced shear and shock-induced tur- 
bulent motions - can not be active in the com- 
parable zone in M51 with high Td ep (§ 15. 1| . Not 
only is shear reduced in this zone, but the gas 
is not characterized by extraordinary line-widths, 
implying that the fewer collapsing clouds in this 
zone are not the result of preferentially enhanced 
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Fig. 8. — (Left) Star formation surface density Hgpji versus observed molecular gas surface density S/f2 measured in 
2.4" radial bins from the center to the edge of the PAWS field, with point colors and shapes as in Figure [4] The size of 
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turbulence. 

Dynamical pressure can also be invoked, as an 
alternative to shear and shocks, to explain the sup- 
pression of star format ion in the Galactic center 
(|Longmore et alj l2012f ). The observed gas kine- 
matics in the Milky Way's central molecular zone 
(CMZ) are consistent with characteristic motions 
along bar orbits, but the evidence for inflow mo- 
tions along a bar shock are less clear; the fall- 
off in the signature 'parallelogram envelope' with 
increasi ng distance from th e center in the pv di- 
agram (jBinnev et al.lll99ll) is not especially ob- 
vious. This could be due to the data quality 
(sensitivity issues), or it may indicate that the 
weak oval nuclear bar does not sustain shear and 
shocksQ In this case, if shear is not present to 
stabilize clouds, dynamical pressure due to non- 
circular motions in the bar may be responsible for 
altering the ability of clouds to form stars. Es- 
timates of the critical surface densities for shear 
and dynamical pressure (see Appendix [D]) could 
help further establish the likelihood of these two 
scenarios. 

5.3. Implications for local and global star 
formation relations 

As demonstrated in Section T4.2.31 the influence 
of dynamical pressure manifests itself by introduc- 
ing departures from a fiducial gas depletion time 
T~dep,o, with higher streaming motions leading to 
longer Tdep- According to the formalism there, this 
fiducial Td ep ,o is representative of virialized clouds, 
and our fitting in the case of M51 (§ I4.2.3[) suggests 
that Tdep,o can be as low as ^0.2Gyr. By analogy 

18 Note that shear is expected to be zero in the absence 
of shocking due to gravitational or viscous torques on the 
gas within the zone of the bar. 



with lKrumholz. Dekel k, McKeel (|2012t l. we view 
Tdep,o as set by the cloud free fall time £// mod- 
ulo some roughly constant dimensionless efficiency 
tff (~0.01 in a supersoni cally turbulent medium; 
iKrumholz fcM cKce 2005[) measuring the intrinsic 
efficiency of star formation in a cloud. But this 
efficiency is modulated by the dynamical pressure 
term so that 



' aep — ° 



(ii) 



This change to the efficiency reflects the influence 
of the environment of a cloud on its star forma- 
tion properties. Alternatively, we can view the 
exponential term as altering the time required to 
compl ete collapse. This would b e mor e similar 
to the lKrumholz. Dekel fc McKee I (|2012t) picture, 
in which the free-fall time depends on the cloud 
'regime'; clouds that are dynamically coupled to 
their environment with a low density contrast have 
shorter depletion times than clouds that are de- 
coupled, at high overdensity. The environmental 
dependence parameterized by this bimodal free- 
fall time leads to a successful fit to observations 
ranging from star forming clouds in the Milky 
Way to high-z starbursts. Below we show that 
our predictions can also explain the observations 
by presenting a form of the star formation law 
that leads to analogo us environment- dependence. 
In co ntrast to the IKrumholz. Dekel fc McKee I 
(|2012l ) picture, our environmental dependence is 
smoothly parameterized by the dynamical pres- 
sure, which affects the coupling to the environ- 
ment by reducing the cloud surface pressure. 

5.3.1. Impact on local trends 

With the change in the depletion time in eq. 
(jlip . the star formation relation can be expressed 
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as 



^SFB, 



^gas (~f + l)vj/4o- 2 



Tdepfi 



(12) 



We see immediately that the range in streaming 
velocities both within and among spiral galaxies 
should contribute to the scatter in th e Kennicutt- 
Schmidt star formation relation (i.e. Leroy ct al. 
2009; Bi gielet all 12001 Fl This is illustrated in 
Figure[5J The left panel shows measurements from 
annuli within M51, while these annuli are plot- 
ted together with integrated measurements from 
Kcnnicutt (1998) for normal spiral galaxies and 
starburst environments on the right. 

Regions in M51 with the highest streaming mo- 
tions sit at the longer gas depletion times, falling 
below the points with relatively less streaming. 
Spiral galaxies with similar non-circular motions 
will show a similar spread in (local) depletion 
time. But more importantly, given that typical 
maximum streaming velocities are 10-15 km s _1 , 
we expect that spiral galaxies will preferentially 
show roughly the same maximum Td ep , especially 
since measurements on the scales of 1 kpc and 
larger tend to be weighted to the bright, high den- 
sity spiral arms. Galaxies with little to no stream- 
ing, on the other hand, will show lower depletion 
times; the gas in these systems will form stars 
closer to the fiducial gas depletion. 

In this light, we might vie w the gas depletion 
time Tde P ^2 Gyr measured bv lBigiel et al.l J2008) 
for a sample of spiral galaxies (ranging in stel- 
lar mass from ^10 7 to 10 11 Mq) as characteris- 
tic of the depletion time in the presence of spi- 
ral streaming motions. At the same time, we can 
attrib ute the lower ~1 Gyr gas depletion times 
in the iKennicuttl (|1998f ) starbursts and also char- 
acteristic of low-mass syste ms assuming a G alac- 
tic conversion factor Xco (|Lerov et al.ll2013t but 



see discussion about the impact of Xco below) 
to the fact that this star formation is occurring 
in a relatively stable dynamical environment. In 
the former case, the star formation typically oc- 
curs at resonances where the streaming motions 
are zero, while in the latter case, streaming is ab- 
sent because these systems fail to develop the non- 
axisymmetric structure that drives these motions 
(the low mass means the disk never becomes un- 
stable). 

5.3.2. Impact on global trends 

In reality, we expect to see a range of gas deple- 
tion times intermediate between these cases. Since 
disks of a certain mass will be unstable to struc- 
ture of a particular wavenumber and pitch angle 
i p , as upheld in the density wave paradigm, we 
can express the i v dependence in terms of galaxy 
mass, or circular velocity V c , using the disper- 
sion relation for m-armed density wave spirals (i.e 

19 Pixel-by-pixel comparisons between S gas and SgFR 
should reveal such departures from a universal gas deple- 
tion time most clearly; the T^ ep measured from larger-area 
averages is the true T^ £p weighted by Hsf R and will thus 
less strongly track low star formation environments. 



Binn ev fc Tremainelll987ft . We use the perturbed 
continuity equation 



(ft — ftp) tanip S 
m Sq 



(13) 



to express the angular velocity difference $7 — Q, p in 
terms of the radial streaming velocities v r , where 
rip is the wave/spiral pattern speed and S (So) is 
the perturbed (unperturbed) gas surface density. 
In this way, we find 
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where V c =flR is the rotational velocity. This re- 
produces the well-known observation that spiral 
arm pitch an gle increases with galaxy mass, V c , or 
Hubble type JKennicuttlll981P i. 

Now from eq. (|13l) we predict that the magni- 
tude of the streaming motions 



v s 



V r (X 



K_S 
m So 
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to lowest order in V c . (Note that this applies 
globally, away from corotation, so that v r ~v^ and 
v s ~v r .) 

Less-tightly wound spirals in the more massive, 
early type disks drive stronger streaming motions. 
This implies that spiral galaxies with higher mass 
(larger V c ) should exhibit longer gas depletion 
times. This is in fact what global measurements of 
Ssf_r and S gas suggest: the normal spiral sample 
in the KS plot on the right in Figure [5] indicates 
an ordering between depletion time and morphol- 
ogy/Hubble type. (Note that while the starbursts 
lie above M51 and the earlier-type spiral galaxies, 
they are more similar to the lower mass later-type 
galaxies, which exhibit similarly short depletion 
times.) This is even clearer in the HERACLES 
sample of galaxies, where integrated measures of 
the gas depletion ti me smoo thly decline from early 
to late spiral types (jLerov et al.ll2013D . 

This dependence of vs on V c accounts for the 
observed weak correlation between T de-o and stel- 
lar mass a mong galaxi es observed by | Lerov et al.l 
(2008, 2009) and Saint onge et"ai] fl2012j). Accord- 
ing to eq. (|13p or (fTS"]) , the scatter in that relation 
at fixed stellar mass can be attributed to, e.g., 
morphologic variations (number and strength of 
arms, pitch angle) and gas content (see below). It 
can also explain the modest success of the 'dynam- 
ical Schmidt' law in fitting observations, like other 
similar prescriptions that express the gas depletion 
time in terms of V c (i.e. lTanll2010l and references 
therein). Here, though, we predict a non- linear 

dependence Tdep ex e c /cr that can be approxi- 
mated by V c only at low V c (and in particular at 
low V c /a). 

5.3.3. High z star formation relation 

Our form of the star formation relation is also 
capable of linking local and high-z star formation. 
Observations suggest that star forming galaxies 
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at earlier times form stars at high e fficiency com- 
pared to their local counterparts (i.e. iDaddi et al.l 
2010), with depletion times comparable to those 
in local starbursts environments. (Local and z—2 
star forming disks fall on a non-lin ear KS rela- 
tion with power-law index n^l.3, e . g. IDaddi et al.1 
[20101: iKrumholz, Dekel fc McKeelfeoia ) High-z 
starbursts and mergers are inferred to consume 
their gas even faster, giving the app earance of a 
bimodality in the K S relation at z^2 (jDaddi et al.l 
2010): IGenzel et aLll2010l ). These authors suggest 
that these two modes of star formation can be 
linked by adopting a form for the star formation 
relation that incorporates the dynamical time, i.e. 
letting the the dyna mical time set the gas deple- 
tion t ime. However, IKrumholz. Dekel fc McKee I 
(2012) demonstrate that this form of the star for- 
mation relation can not be consistently applied 
to local systems. Instead, the fiducial gas deple- 
tion time is argued to be lower in high-z starbursts 
likely due to the expected high densities that re- 
sult under the compressive weight of the ISM. As 
is clearer below, our picture presents a parameter- 
ization of this weight in terms of V c /a or the gas 
mass fraction f g . 

To more clearly illustrate the differences we ex- 
pect to see between systems at low and high z, 
we alternatively express the streaming velocity in 
terms of the Toomre Q and gas-to-total mass ratio 
f g . Using that 



Q = 



y/2(P+l)an 



(16) 



where /3=d\nil/d\nR—l for solid body rotation 
and j3=0 in the case of fiat rotation curves, 
V^^TrCStotRedge, and assuming that the stellar 
and gaseous disks cover roughly the same area so 
that fg=T, gas /T, to t then, together with eq. (|T5|) . 
we find 

S 1 v/2(/3 + l) 



v s 

a 2jq rn 



f 9 Q 
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Star forming disks at high-z have lower gas de- 
pletion times than similarly marginally unstable 
disks at z—0 due to the higher gas fraction at 
earlier times. High-z starbursts and merging sys- 
tems can have even lower depletion times, given 
the likely higher gas fractions. 

Note that the observed trends also follow from 
eq. (I15p . as the ratio V c /a is lower in star forming 
disks at earlier times, and even lower in starburst- 
ing mergers (i.e. LIRGS/ULIRGS), given the en- 
hanced turbulence and larger a. This may explain 
why a linear dependence of 'Ssfr on V c can be 
used to fit observations at high z, where the quan- 
tity Vc/cr is intrinsically low, and thus a suitable 
approximation of the exponential dependence we 
find (see § !5.3.2p . We emphasize, though, that V c 
is only part of the description, as the ratio V c /a 
provides the more fundamental measure of the im- 
portance of dynamical pressure. As described in 
§ 14.11 even when streaming motions are present, 
dynamical pressure is not important if clouds can 



equalize before they undergo translation by more 
than a cloud length, that is, if a>vs- With eq. 
(I15p . this is equivalent to the condition E 9as ~S toi 
for gas disks in hydrostatic equilibrium (recall that 
by writing eq. (|15[) . we replace the existence of 
streaming motions by the requirement for a disk 
massive enough for instabilities to form). So even 
in the presence of streaming motions, a high gas 
fraction prevents a lengthening in the depletion 
time. 

Consider, for example, the difference between 
starbursts at low and high z, which each have 
short depletion times but very different V c . In 
the case of local starbursts, the high gas densities 
(and high star formation) are possible because of 
the absence of motion relative to the background 
rotating potential; locally motionless gas (fed from 
neighboring regions) can retain and build a molec- 
ular reservoir. This scenario is different from the 
conditions that lead to bursts of star formation at 
higher redshift, where accreted external gas builds 
up globally high surface densities and streaming 
motions themselves ca n be quite high (i.e . given 
the high measured V c ; IGenzel et al.ll20l6T ). Even 
though these streaming motions are present, dy- 
namical pressure is less important, according to 
the measured ratio V c /a. 

5.3.4. Comparison with other findings 

In our picture, shorter depletion times are asso- 
ciated with greater levels of turbulence and higher 
a. This is also a feat ure of the model proposed by 
iRenaud et al.l (|2012f ) based on feedback-regulated 
turbulence. The higher turbulence in z—2 merg- 
ers as compared to z—2 disks leads to a higher 
fraction of (locally) dense gas, higher overall star 
formation rate and thus a vertical offset in T,sfr 
vs. Tigas space. Despite consistency with this de- 
pendence of Tdep on <t, tho ugh, we interp r et the 
role of a differently. In the IRenaud et aLl {2012) 
model, lower levels of turbulence lead to less ef- 
fective turbulent-triggering of star formation and 
longer depletion times. In our picture, a lower a 
leads to a more destabilized system and thus to the 
onset of streaming motions, which lengthen the 
depletion time. Note that both of these pictures 
appear to imply that star formation consumes gas 
quicker when the gas is more stabilized (accord- 
ing to the increase in Toomre Q with the higher 
effective gas sound speed in the presence of greater 
turbulence), perhaps contrary to expectations. 

Interestingly, at fixed a, eqs. (fTTj) and (fTT]) 
together predict that the gas depletion time in- 
creases nonlinearly with decreasing gas fraction, 
or that the star formation efficiency increases non- 
linearly with gas fraction. In addition, more cen- 
trally concentrated gas, situated in the portion of 
the disk where the rotation curve is rising and 
/3~1, will have slightly longer Td ep than more 
extended disks with /3=0. In the former case, 
shear prevents gas from forming collapsing struc- 
tures. Our pic ture is thus also able t o describe 
the findings of iSaintonge et al.l ((2012T ) and may 
prove descriptive, from an observational perspec- 
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five, of 'morphologi cal quenching' as envisioned by 
iMartig et all (|2009 ). 

We can now put this all together to predict a 
dependence 

log EsiTfl = log £ 9as - log t///e// — — (18) 

for a Toomre unstable disk (Qwl), and taking for 
simplicity Ti/YiQ^m—2 an d 7=-1.5 for the cloud 
mass spectrum index (i.e. Rosolowsky 2005). As- 
suming a universal tff/e/f that sets the intercept 
of a 'fiducial' linear relation between log T,sfr and 
E ffQS , then systems will fall progressively below 
this line the smaller their gas fraction. This form 
of the star formation relation can consistently de- 
scribe the strong offsets between various popula- 
tions in the standard KS plot as assemb l ed an d 
presented by iKrumholz. Dekel fc McKee I (|2012t ). 
Systems with the highest gas fractions - high-z 
starbursts, local starbursts, and even including 
star-forming clouds in the Milky Way - are off- 
set to the shortest depletion times. Below these 
fall the high-z star forming disks at intermediate 
Tdepi followed by the local star forming galaxies at 
the longest Tdep- 

Of course, streaming motions may not be re- 
sponsible for the entirety of the offset in observed 
gas depletion time both within and among galax- 
ies. Some part of the offset can be attributed to, 
and removed by, variation in the CO-to-H2 conver- 
sion factor (Xco)- The bimodality between local 
and high-z systems is arguably enhanced due to 
the a doption of a bimodal Xco ([Narayanan et all 
120 111 ). Allowing for a dependence of Xco on 
metallicity, low- mass dwarfs (with high Xco) ap- 
pear to have s imilar gas depletion times as more 
massive disks (|Schruba et alj|2012t unlike in Fig- 
ure [8] where a single Xco is assumed) . 

On the other hand, it seems that not all de- 
pletion time variations between different popula- 
tions can be minimized in this way. In M33, a 
spiral galaxy intermediate in mass and metallic- 
ity between the MW and Local Group dwarfs, 
the depletion t i me is shorter than for spirals in 
the IKennicuttl (11998) sample, even adopting a 
slightly higher t han Galactic Xco appropriate for 
its metal licity (iBraine et al.ll2011[ ). In addition, 
although iLerov et al.l (12013) find that most varia- 
tion in Tdep internal to galaxies can be accounted 
for by varying Xco, the differences in Xco be- 
tween the di sks and central regions of galaxies 
measured by ISandstrom et al.1 (J2013I ) in fact un- 
derscore the ex istence of two diff erent 'modes' of 
star formation ()Lerov et alj|2013l ). 

6. CONCLUDING REMARKS 

In this paper, we examined the influence of 
non-axisymmetric stellar structure (i.e. nuclear 
bar, spiral arms) on gas flows, molecular cloud 
properties and star formation in the inner disk 
of the iconic spiral galaxy M51. Leveraging 
our unique view of gas motions, which includes 
both measurement of present-day torques and 



non-circular streaming motions decomposed from 
the line-of-sight velocity field, we establish an 
anti-correlation between gas flows and strong star 
formation in M51. To explain the observed gas 
flow and star formation patterns we developed a 
simple model, in which dynamical pressure is a 
critical parameter for determining cloud stability 
against gravitational collapse. We report the 
following results and conclusions: 

1. A radial profile of the stellar torques across 
the PAWS field reveals distinct dynamical zones 
within which gas is driven radially inwards or 
outwards in response to the non-axisymmetric 
structure present in M51's inner disk. In the nu- 
clear bar region (< 20") and inside the corotation 
resonance of the first spiral arm (35" < R < 55") 
gas is driven strongly radially inwards. Elsewhere 
in the disk gas is either stationary or experiences 
a radially outward torque. 

2. We use a unique method to decompose the 
line-of-sight velocity field of the CO emission 
into its radial and azimuthal components. Across 
the PAWS field, the sign and magnitude of 
our derived radial gas streaming motions agree 
very well with our expectations for how the 
molecular ISM should respond to torquing by the 
non-axisymmetric stellar structure. 

3. Comparison of the radial torque profile and 
gas streaming motions with the molecular gas 
surface density and star formation rate (as traced 
by the combination of Ha and 24 fim emission), 
reveals a clear anti-correlation between strong 
gas flows and active high-mass star formation. 
In M51, radially inflowing gas appears to be less 
efficiently forming stars. More generally, gas 
sitting near the corotation of the disk and the 
bar or spiral structure forms stars more efficiently 
than gas that is far from corotation and in motion 
relative to the background potential. 

4. We propose that the complementary pat- 
terns of gas flow and star formation in the 
PAWS field are evidence for the importance of 
external pressure on the properties and stability 
of star- forming clouds in M51. More precisely, we 
develop a simple model based on the Bernoulli 
principle that describes changes in cloud surface 
pressure in the presence of strong streaming 
motions (dynamical pressure). Such changes in 
surface pressure can make the difference between 
stability and collapse when, as in M51s spiral 
arms, clouds are dynamically coupled to their 
environment (i.e. when the cloud internal and 
external pressures are similar). In M51, cloud 
stabilization occurs preferentially as the result of 
strong radial inflow motions in regions of negative 
torque, but we expect that outflow motions may 
serve the same role in other systems. 

5. Our model successfully reproduces sev- 
eral key observations. We can retrieve the 
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environment-dependent variation of the slope 
of the cloud mass spectrum measured directly 
from the cloud distribution. We can also explain 
the overall difference in the masses of clouds 
between the arm and interarm. Sensitivity to 
changes in dynamical pressure are also consis- 
tent with the observation that not all of M51s 
clouds are virialized, as evidenced by the lack 
of a clear size-linewidth relation for clouds in M51. 

6. We investigate other potential sources of 
cloud stability (shear, galactic shocks, stellar 
feedback-driven turbulence) and find that they 
can not uniformly explain the observed non- 
monotonic radial dependence of the gas depletion 
time. Cloud formation mechanism also has little 
influence on the ability of clouds in M51 to 
collapse and form stars: while cloud collision 
and agglomeration in the spiral shock are likely 
present, we find no obvious relation to star 
formation in M51s current gas reservoir. 

7. Although non-circular motions in M51 are 
particularly high, the influence of dynamical 
pressure should be common-place in spiral galax- 
ies. We suggest that gas flows within galaxies 
are a source of scatter in the Kennicutt-Schmidt 
star formation relation both within and among 
galaxies. Late- type spiral disks with tightly 
wound spirals that drive weak streaming motions 
will form their stars at a faster rate per unit 
gas mass than early-type disks where streaming 



motions are higher. 

8. We propose that our model for the depen- 
dence of gas depletion time on dynamical pressure 
links low- and high-z star formation from dense, 
molecular material. In particular, we suggest 
that systems form stars with progressively shorter 
depletion times the higher their gas fraction. 
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APPENDIX 
A. MASS-BASED ROTATION CURVE FOR M51 

We use the standard approach to assemble an estimate of the rotation curve in M51 by summing the 
individual contributions of the stellar and gaseous components to the circular velocity, adopting a simple 
mass-follows-light model appropriate for spiral galaxies (i.e. the rotation curve can be accounted for with 
baryons, alone) PI 

This requires first constructing radial profiles of the stellar and gas surface densities. The total gas 
surface densit y T, nns is assem bled from measurements of Y^hi and Sg „ traced by the T HINGS HI zeroth- 
moment map (jWalter et al.l f2008) and PAWS CO(l-O) observations (|Petv et al1l2013| ). respectively, and 
includes a factor of 1.36 to account for the presence of Helium. As in §[2] we assume a constant CO-to-i?2 
conversion factor of Aco=2x 10 20 cm -2 (K km s -1 ) -1 and also include a sky-to-disk plane correction of 
cos 21°. 

The stellar surface density £* is traced by the S 4 G 3.6 fim map of the old stellar light (corrected for 
non-stellar emission; Meidt et al. 2012). The (assumed) global 3.6 fim mass-to-light ratio T 3 . 6 (with 
which E Jr =T3.6i3.6COs21°) is the single free parameter to be fit in comparing the modeled and observed 
circular velocities. The final surface density profiles are shown in Figure [9] Below we adopt models of 
the individual exponential disk and Gaussian bulge components of the stellar surface brightness. 

To estimate the circular velocities associated with the gaseous and stellar disks we adopt the thin disk 
approximation and model the bulge separately. The circular velocity at radius R (in the plane of the 
disk) is given by 



V? D (R) = 2irGR 



Z(R')J (kR')J 1 {kR)R'kdkdR' 



(Al) 



where Jo and J\ are cylindrical Bessel functions (e.g. iBarnes fc Sellwoo"c3 2003; B arnes fe S ellwood 2004). 
Without accounting for disk thickness, which would reduce the circular velocity, the derived T is slightly 

underestima ted. 

Following iPalunas fc Williams! (|2000l ) we estimate the contribution from the stellar bulge to V c by 

20 The maximum disk hypothesis h as been foun d to yield good fits inside the optical radius of 3/4 of the galaxies in 
spiral sample studied by Palunas & Williams (2000). 
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Fig. 9. — Surface density profiles in the central 9 kpc of M51. The thin gray dotted line shows the total gas surface density 
-T,ff 2 +T,j{j. The solid gray and black lines show the modeled and observed stellar surface density E* assuming the 



best-fit stellar M*/L3.6 (determined below), 
(dashed) and Gaussian bulge (dot-dashed). 



The model components are shown in thin gray lines: exponential stellar disk 



modeling the bulge component of the stellar surface brightness profile with a Multi-Gaussian expansion 



HB{a B ) = \/~£fc =1 — 
7T R B 



^ e ~al/R\ 



(A2) 



where clb=R 2 + (z/(l — £b)) 2 is the bulge semi-major axis length, eB = l-q P ,B is the bulge ellipticity, q Pi s 
is the apparent bulge axis ratio, f=cosi 2 + smi 2 /q B and, for each of k components, Ck is the total light 
and R Bk is the scale length. 

This 2D elliptical Gaussian distribution deprojects to a 3D spheroidal distribution with intrinsic axis 
ratio qB—iqi b ~ cos « 2 )/ sin i 2 whereby (in the plane of the disk) 

rR v(R' 2 )R' 2 dR' 



V' B (R) = 4irG?a. 6 Jl 



B 



(I 



VR 2 



e|i?' 2 



(A3) 



with &b the bulge eccentricity, which we take to be zero for simplicity (given the small degree of asymmetry 
in the central projected light distribution) . We expect this assumption to have less consequence than our 
requirement for a global T3.g, uniform from bulge to disk. 

Modeled circular velocities are plotted in Figure [TU] together with two independent esti mates of the 
rotational velocity measured from the PAWS CO (1-0) observations (|Colombo et al.ir2013al) . The thick 
black curve shows the rotational velocity measured by fitting tilted-rings to th e line-of-sight vel ocities 
in the combined interferometric and single dish data (with native 1" resolution; iPetv et a l. 2013) using 
the GIPSY task ROTCUR with fixed center, position angle PA=173° and inclination angle i=21°. The 
thick red curve shows the rotational velocity measured from the 30m single-dish data (12" resolution) 
assuming the same kinematic parameters. We expect this curve to be a fairer representation of the 
true circular velocity, as beam-smearing at lower resolution minimizes the contribution from non-circular 
streaming motions to the measured rotational velocity (evident as wiggles in the thin curve from the higher 
resolution data). Our optimal T3.6 is chosen based on agreement with this estimate of t he circular velocity, 
as de scribed below. We adopt a smooth three-parameter fit to the total velocity (i.e. iFaber fc Gal lagher 
1979) as our final model for V c . With our mass-follows-light model w e find a best-fit stellar T3 «=0.4 5 
±0.15. This value agrees well with the value T 3 . 6 =0.5 measured by lEskew. Zaritskv fc Mei dt (2012), 
which is its elf very near the y alue that has been previously assumed to convert between 3.6 /im luminosity 
and mass ([Lerov et al. 200cf). The maximum disk hypothesis is evidently reliable in this case, at least 
throughout the PAWS field of view. 
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Fig. 10. — Rotation curves in the center (ij<4.5 kpc) of M51. Circular velocities derived from the PAWS CO(1-0) data 
are shown in thin lines (combined data, black; single-dish, red; errors are omitted for clarity). The thick solid black line 
shows the total V c , calculated as the sum of the individual modeled circular velocities shown in thin gray lines: stellar bulge 
component (dot-dashed), stellar disk component (dashed), gas (dotted). Comparing the red curve with our smooth fit to 
the total circular velocity (thick gray dotted line) gives a best-fit M/L=0.45±0.15. 




Fig. 11. — Map of torques in the inertial frame R X V<J>, generated from the 3.6 fim map of the old stellar light tracing the 
stellar potential $ (see text). White (black) corresponds to positive (negative) torques that drive motions radially outward 
(inward). Contours of the CO intensity are overlayed in red. The green bar at right indicates 40". Image dimensions: 
~300"x200". 

B. DEFINING THE DYNAMICAL ENVIRONMENT: PRESENT-DAY TORQUES 

The impact of environment is underlined by present-day torques exerted by non-axisymmetric struc- 
ture in the disk. We use information supplied by these gravitational torques to define the dynamical 
environment. 

Gravitational torque estimation 

The S 4 G 3.6 \xm map of the old stellar light (Meidt et al. 2012; from which a contribution from 
non-stellar emission has been removed) presents us with an optimal view of the backbone of present-day 
gravitational torqueso To generate a map of the stellar potential we first deproj ect the image according 
to the inclination i—21° and major axis position angle PA=173° determined bv lColombo et all (|2013bf ) 
from gas kinematics. We then apply a constant stellar mass-to-light ratio (M/L) to the 3.6 fj,m image 



21 Note that the influence of the companion is included to the extent that it has contributed to the reorganization of the 
stellar mass and structure in the disk mapped in the stellar potential. 
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and map this into density by assuming a simple tanh vertical distribution with uniform scale height. Like 
the stellar M/L, the adopted vertical distribution will affect the magnitude of the measured torques but 
not the radial dependence highlighted below, which is our primary concern. Triaxiality in the DM halo 
could potentially introduce a change in the radial behavior, but we expect this to be minimal over our 
relatively small and centralized field of view covering the inner 9 kpc. The dominance of the stellar mass 
over this area is confirmed by our study of gas kinematics in § 13.21 which presents supporting evidence 
for the conclusions drawn from the stellar potential alone. 

From the potential we calculate a map of the quantity -R x V<i> (Figure fTTI) represent i ng th e in- 
ertial torque per unit gas mass, following iGarcia-Burillo et al.l (J2005D . IGarcia-Burillo et all (J2009D and 
lHaan et al.l ([2009). In the rotating (non-inertial) frame forcing includes a contribution from the 'fictitious' 
centrifugal and coriolis terms. Since the cross product of radius R with these terms go identically to zero 
in the former case and to zero with averaging over to 2n in azimuth in the latter, then a measure of 
the time average of the rate of change of the angular momentum, the torque, can be obtained by taking 
azimuthal averages of R x V<1>, with we i ghting by t he molecula r gas s urface density. This is equivalent 
to the approach of IGarcia-Burillo et all (|2005l) and lHaan et all (|2009fl where the molecular gas surface 
density represents the probability of finding the gas at its current position at present. Here we use the 
PAWS observations of CO(l-O) emission to trace the molecular gas surface density assuming a CO-to-Ff2 
conversion factor Aco=2xl0 20 cm~ 2 (K km s -1 ) -1 (|Colombo et al.ll2013al ). These observations include 
single dish data, which has been found to be a critic al addition to the gas tracers used for this purpose 
(|Garcia-Burillo et al.l 120051 : Ivan der Laan et all 120111 ) ; extended emission on the largest scales must be 
included for accurate determination of the net torque. Note that the assumed conversion factor is another 
potential influence on the magnitude of the derived torques. For this reason in the upcoming section we 
present the torque profile in units of the net torque over the PAWS field of view. 



Inventory of distinct dynamical environments 

Figure [T] in the main text shows the radial profile of azimuthally averaged torques across the PAWS field 
of view. Each crossing from negative to positive marks a corotation radius (CR), while crossings from 
positive to negative torque coincide with a switch in predominance to a new pattern. Three dynamically 
distinct zones are identified, as described below. 

The first CR in Figure Q] occurs just outside the " butterfly" pattern in the c entral 20" of the torque 
map, which is characteristic of the bar influence (i.e IGarcia-Burillo et al.1 119 93) . Comparison with the 
velocity field confirms that the gas is responding to torquing by the stellar bar (first identified in the 
NIR by iZaritskv. Rix "fc Rickc 1993), even though the stellar bar lacks a molecular counterpart. The 
bar is responsible for the +10-15° twist in the orientation of the line-of-nodes, previ ously interpreted a s 
evidence for variation in the galaxy major axis tracing a warp or twist of the disk (jShettv et al.l 120071) . 
This marks the first evidence that the bar, rather than the spiral, dominates the kinematics in the central 
30". 

The influence of the bar also extends to gas morphology, as evidenced by the build-up of a molecular 
ring at the overlap with the molecular spiral arms just exterior, near R=35". Opposing torques from 
inside the ring, where the bar drives gas outwards, and outside the ring, where torques drive gas inwards 
(see below) , result in the pile- up of gas at the location of the ring. 

Beyond the second corotation at Rcr— 55" gas is again driven outwards, across the outer zone of the 
CO-bright spiral arms. Radial outflow continues until the start of the third distinct pattern, near R=85", 
where the direction of flow is again reversed. This loc ation is consistent with previous estimates for the 
start of the distinct outer, material spiral pattern (e.g. lTulry||1974t iVogel et alj|1993t iMeidt et ai]|2008bD 
which continues beyond the edge of our field of view. 

In total, the gravitational torques expose three dynamically distinct zones dominated by at least three 
unique non-axisymmctric structures. The number, multiplicity and r adial domains of patte rns present 
in the disk will be the subject of a more detailed study in the future ([Colombo et al.ll2013bl ). Below we 
briefly review supporting evidence for our identification of these three main environments. 



Complimentary evidence for multiple patterns 

We expect at least three distinct pattern speeds in the central 9 kpc of M51, coinciding with three 
distinct structures: the bar, the main spiral and an outer spiral. For the main spiral pattern with 
corotation at Rcr=55" we estimate a pattern speed of £l p ~90 km s _1 kpc -1 based on inter section with 
our ado pted angular frequency curve (see Fig. [T^]). T his is very close to the value measured by IMeidt et all 
(|2008bl) with the radial Tremaine- Weinberg fTWR: iMerrifield. Rand fc Meidtll200l IMeidt et aH l2008a) 
method using lower-resolution molecular gas observations as a kinematic tracer, and assuming a similar 
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Fig. 12. — Pattern speed estimates for structure in the disk of M51. Angular frequency curves are shown: Q (Black), 
SI — re/4 (Gray), Q ± re/2 (Light Gray). Estimates for the bar and spiral pattern speeds implied by the corotation radii 
revealed by the torque profile in Figure [T] are shown as gray horizontal lines. A possible second spiral speed (see text for 
discussion) and the rotation of the material arms are also indicated in gray horizontal dashed lines. 

major axis position angle of PA=175°E3S 

At smaller radii, the gravitational torques suggests a transition to a distin ct bar pattern, with Rcr— 20", 
which is comparable to the location recently suggested bv lZhang & Butal (|2012[ ). From this we estimate 
a bar pattern speed of $7b~200 km s -1 kpc -1 . In this case, the bar and spiral pattern speeds offer a 
physically plausible scenario, as the corotation of the bar coincides with the inner 4:1 r esonance of the spi- 
ral. Resonance overlaps of this kind have been ide ntified in barred spiral simulations (Rautiai nen fc SaTol 
1999; Debatti sta et al.l|2006t JMinchev et al.ll2011f ) and may be characteristic of nonlinear mode coupling 
(e.g.. rTagger et alJ ll987 and Sv gnet et al.lll9990 . whereby energy and angular momentum are transferred 
between the modes. Note that with this set of pattern speeds, the molecular ring sits near the spiral's 
inner 4:1 resonance (and near the bar's outer 4:1 resonance) and forms at the intersection of positive 
torquing by the bar outside its CR with the negative torquing by the spiral inside its CR. The inner 
4:1 resonance has p reviously been suggested to favor the population of gas in the form of such a ring 
(Garcia-Burillo et al. 2009). (We also note that the coincidence of the end of the bar with its CR makes 
this bar a so-called 'fast' bar.) 

This same process could be at work further out in the disk, in the zone between 35" and 70". At its 
Rcr=55" the speed J7 P =90 kms^^^kpc" 1 could overlap with the inner 4:1 resonance of a lower speed of 
^p,2~55 kms~ 1 kpc~~ 1 (see Figure IT2l . just before the start of the material arms. This second spiral speed 
falls just within the range of speeds suggested for a similar second speed in the TWR solution at 170° 
measured by Meidt et al. (2008). 

The existence of this second more slowly rotating pattern may be plausible for several reasons: 1) the gas 
kinematics less clearly show the pattern of outflow in the post-corotation zone 55"<i?<70" as would be 
expected for a pattern outside its CR0 Two patterns, each propagating inside their CR would maintain 
radial inflow through this zone. 2) The spatial offset between young star and gas tracers does not clearly 
switch sign across the corotation at Rcr— 55" as predicted in response to a spiral density wave around 
CR. Again, two patterns, both propagating inside CR, and both driving gas radially inwards, removes 
the possibility that young stars will be found upstream of the gaseous spiral, as otherwise expected 
for a single pattern outside its CR. 3) While the 3.6 p,m map of old stellar light shows a pronounced 

22 Accounting for differences in adopted distance and inclination angle 

23 Two speeds from an alternatively plausible TWR solution, which assumes a lower PA=170°, straddle our adopted 
v alue Sl p ~90 km s -1 kpc -1 . Based on evidence here for the bar influence on m orphology and molec ular gas kinematics (see 
£j 13. 21 1. though, it seems likely that the higher inner pattern speed measured by Meidt et al. (2008b) reflects a combination 
of the speed a t n p =90 km s k pc - and a yet higher bar pattern speed. This circumstance has been encountered in 
other galaxies (Meidt et al. 2009), where the resolution in the TWR solution was limited by that of the kinematic tracer. 
According to our estimate of the bar CR, the bar zone would have been probed by only the inner three to four ~6" radial 
bins used by Meidt et al. (2008b). In addition, given the orientation of twisted line-of-nodes manifest by the bar, the higher 
PA=175° major axis orientation for the TWR solution presumably disfavors the detection of the bar influence, more so than 
for a lo wer assumed PA; as for the TWR solution at PA=170° a third solution with PA=165° measured by Me idt et al.l 
(2008b) suggests an inner pattern speed >90 km s -1 kpc -1 . 

24 The radial torque profile <T>(R) (Figurc[T} suggests that radial gas inflow stops at Rcr=55" , rather than continue 
until R=70" , as suggested by the radial streaming motions v r , measured in Appendix [C] This difference could be alleviated 
with a constant offset to either v r or <T> (with an impact on the locations of the other CR radii, as well), or it could 
indicate a more complicated gaseous response to the torquing. 
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Fig. 13. — (Left) CO intensity-weighted velocities parallel (vn, blue) and perpendicular (v±, red) to the spiral arm, 
calculated in radial bins as described in the text from the reconstructed radial and azimuthal velocities (see Figure [2] in 
the main text). Dashed lines represent the uncertainty associated with each velocity component, calculated from the errors 
shown in Figure [2] (i.e. the rms dispersion in solutions with ±5° variation in the assumed spiral arm pitch angle i p and the 
uncertainty introduced by our adopted rotation curve). The intensity- weighting focuses our measurements on the CO-bright 
spiral arms. (Right) The difference in the magnitude of non-circular streaming motions in the arm and interarm. 



m=2 perturbation, this is identifiable as a two-armed spiral only outside R^55" . Inside this radius, in 
the zone 35 // <i?<55" , ellipse fits to 3.6 fim isophotes reveal a more oval, bar-like structure (see also 
ISchinnerer et al.ll2013D . Other stellar traces (from FUV to H a ) reveal a similarly marked absence of a 
clear spiral pattern in this zone. In this case, the higher pattern speed of f2 p ~90 kms~ 1 kpc~ 1 might apply 
to a bar terminating at its corotation radius, while the lower pattern speed 17 P! 2~55 kms~ 1 kpc~ 1 could 
describe the strong spiral pattern. We might then interpret the gaseous response to the oval perturbation, 
which is very definitely in the form of a two-armed spiral, as equivalent to dust lanes running along the 
length of a stellar bar. Alternatively, the main spiral pattern could overlap with a lower, outer pattern 
with speed $7 P =20 kms~ 1 kpc~ 1 near where its OLR and the outer pattern's ILR coincide (see Figure 

USD. 

For simplicity in what follows, we will describe the dynamical environment by three distinct zones 
dominated by an inner, or nuclear, bar and two (inner and outer) molecular spirals. As explored in § 
13. 2\ the observed gas motions reflect the predicted combination of radial inflow and outflow that, in the 
spiral arm region, correspond to flow along and through the arm. 

C. DECOMPOSING THE PAWS LINE-OF-SIGHT VELOCITY FIELD 

The nearly logarithmic spiral in M51, spanning over 5 kpc in radius, presents a unique opportunity to 
decompose the line-of-sight velocity 

V[ os (R,9) = V sys 

+ [(V c (R)+v^(R,9))sm(9) 
+ v r (R, 9) cos (9)] sin i 

(CI) 

in to its radial and azimuthal components v r and v$. Here V sys is galaxy system ic velocity, V r is the dis k 
circular velocity, and 9 is measured with respect to the line of nodes. Following iColombo et al.1 (|2C)I3b), 
we adopt an inclination of i=21° and a major axis position angle of 0^=173°. 

The key to decomposing Vi os is the idealized arm frame in which streaming motions v r and v$ will 
be nearly constant in small segments at constant arm phase ip parallel to each arm. Adopting a spiral 
arm pitch angle of i p =21° we transform from cartesian to log-polar coordinates and define a grid of 
rectangular segments or bins aligned with the spiral arm. In each bin we solve for v r and v^ via simple 
matrix inversion, according to 



{Vi os - V sys )/ sini 
-V c sin{6)-- 



■v r (ip) cos(9) 
-v^(t/j)sm(6), 



where the coefficients are specified by the geometric factors. 
We calculate maps assuming three different pitch angles: a nominal i 



(C2) 
,,-21° (i.e. IShettv et alll200l 
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and a second set offset by ±5°. The average of the solutions from these three cases define our estimate 
for each component, while the dispersion in these values supplies a measure of our uncertainty, together 
with the uncertainty introduced by our adopted rotation curve (see Appendix [X]). To accommodate 
for real variation in the molecular spiral arm pitch angle (determined to be 5° between the arms and 
as much as 20° along the arms by Patrikec v et al.l [2006} ) . we size the spiral segments over which the 
solution is calculated to 5.6° wide and Mogr(")=0.011. These segments cover a smaller area toward the 
center, so the solution in the bar zone, where a logarithmic-polar deprojection may be less fitting, should 
nevertheless be accurate. 

Prior to solution the contribution of the circular velocity to the observed l.o.s. velocity field is removed 
(see eq. IC2[) . We adopt the circular velocity estimated from our ID mass model of M51, constructed 
as described in Appendix 1. By using a mass-based reconstruction of the circular velocity we avoid the 
uncertainties typically inherited by circular velocity estimates made directly from the l.o.s. velocity field 
(e.g., with tilted-ring fits), namely the contribution of azimuthal streaming motions. The result is a 
much more smooth, physically motivated model of the circular velocity curve that should insure accurate 
streaming motion solutions. We use the uncertainty introduced in the rotation curve by our choice 
of stellar M*/L3.6=0.45±0.15 (Appendix 1) to define the error in our measured streaming velocities, 
together with the measured pitch angle uncertainty. Both are typically <25%. 

The tangential and radial streaming motions (shown in Figure [2] in the main text) translate into 
motion directed along and through the arm, channeled by the stellar spiral. This is depicted more clearly 
in the top left of Figure [TBI showing intensity- weighted averages of the velocities in the spiral arm frame, 
assuming i p =21°. Here 

V± = v r cos(i p ) + v,p sin(ip) 

I'll = — v r sin(ip) -(- v^ cos(«p) (C3) 

where v± runs perpendicular to the arm and vu runs parallel to the arm. Note that in the zone 40"<r<60", 
the velocities v±_ perpendicular to the spiral arm arise almost completely with the radial motions since 
the azimuthal streaming in this zone is at a minimum. The spiral arm pitch angle is meanwhile increasing 
here, by up to 30°, so that the radial streaming motions could be distributed slightly more in the direction 
parallel to the spiral. This would lead to larger v» than measured here (assuming i p —21°), although the 
perpendicular motions still dominate in this case. 

As described in § 13.21 we use CO intensity as a weighting factor in binned azimuthal averages to reveal 
the velocities characteristic of the (CO bright) spiral arms. Profiles with inverse-intensity weighting, 
which highlight the inter-arm streaming velocities, are qualitatively similar but show velocities ~ 5-10 
km s _1 lower than the profiles in Figure[2] The difference in the magnitudes of streaming motions in the 
arm and interarm is shown in the right panel of Figure 1131 

D. GRAVITATIONAL DISK STABILITY 

In this section we examine in greater detail the stability of the molecule-dominated gas disk in M51. 
We quantify and compare the critical surface density for dynamical pressure to that in the presence of 
shear, Coriolis forces and tidal forces, three sources of cloud stability of potential consequence for the 
global patte rn of sta r form ation- 
Following iToomra (|1964f l we take 

y -— rnn 

7T(jr 

as defining the threshold for stability in a thin rotating gas disk, whereby the Toomre Q—\ when St oom re = 
Ti gas . Here a is the gas velocity dispersion, k is the epicyclic frequency measured from our rotation curve 
model (see Appendix [AJ and a is a dimensionless factor calibrated empiricallyo (I n the following we 
let a = 1 but also discuss the impact of variation in this value.) Below this threshold, Coriolis forces 
spin up the gas so that centrifugal forces counter the self-gravity of the gas, resulting in a suppression of 
cloud (and star) formation. Above the threshold, gas is unstable to large-scale collapse, leading to star 
formation. 

The threshold for cloud stability in the presence of shear, as derived by lElmegreenl ([1993) (see 
iLuna et al.11200 6). can be written 

^shear = 7^ \^^) 

7T(_r 

where A is the Oort parameter and again with a dimensionless factor a a, which we set to unity. We 
consider separately the background shear due to differential rotation and the spiral arm shear, accounting 
for gradients in the azimuthal spiral arm streaming in the definition of A (in addition to background 
shear), as in ? l5.1.1l in the main text. 

25 Values a 7^1 are chosen in orde r that Q=l at the location of the observed star formation threshold. As suggested by 
ISchavei H2004T) and dc Blok & Walter (2006) deviations in a from unity may be commonly required only because measure- 
ments of the gas velocity dispersion from observations of the atomic phase systematically overestim ate the true a of gas in a 
cold phase, which we can better approximate with the PAWS observations (as confirmed in M51 by Colombo ct al. 2013b). 



2(> 



Alternatively, cloud shearing via tidal forces will become ineffective above the threshold 
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following iKennev et al.l (|1993l ). with Oort B = O — A. Note that these formulations of cloud stability 
ignore any modulation/dissipation of shear, tidal or Coriolis forces due to additional effects, such as 
magnetic fields or viscosity. 
According to eq. (5), the critical surface density in the presence of dynamical pressure can be written 
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letting the characteristic surface density of virialized clouds T, V i r =a 2 /ttGH, and using that the gas scale 
height H is approximately c/f2 near to hydrostatic equilibrium in self-gravitating gas. 

In Figure [14] we compare the radial profiles of this set of critical surface densities with the observed arm 
and inter-arm gas surface densities. As shown on the left, the shear critical surface density T, s h ea r falls 
below the observed surface density in the arm almost everywhere, suggesting that shear is not responsible 
for destroying clouds in the arm, consistent with the argument of, e.g., Elmegreen et al. (2003). The 
reverse appears to be true in the interarm where, as considered bv lColombo et ail (|2013a| ). the reduced 
cloud number density (compared to the arm) reflects an efficient cloud destruction mechanism. 

In the zone of the bar, the shear and Toomre critical surface densities approach and even exceed the 
observed surface density, lending support to the idea that Coriolis forces plus shear (and shocks) in 
the bar can stabilize clouds and prevent subsequent star formation. However, several lines of evidence 
suggest that cloud stabilization through this more traditional shear/shock scenario in the bar is weak 
and likely ineffective: 1) the bar itself is weak, as indicated by its small size and oval shape, 2) there is 
little evidence for a clearly defined pair of dust lanes running along the length of the bar (traced by either 
molecular gas or optical extinction), which are characteristic of a strong bar shock and 3) as shown on the 
right of Figure Qj] the threshold for stabilization through dynamical pressure is higher than these other 
thresholds, exceeding T, gas by as much a s an order of magnitude. Indeed, the critical surface density for 
dynamical pressure as expressed in eq. (|D4|) will always exceed the traditional Toomre critical surface 
density where i>s>2o-\/ln 5 for 5 = k/VL (i.e. 8 — v2,l or 2 for rising, flat or Keplerian rotation curves, 
respectively.) 

According to Figure [141 it appears that star formation in the bar zone may be suppressed and related 
most directly to stabilization through dynamical pressure. Stabilization via dynamical pressure is also 
evidently the strongest (if not the only) interpretation for the suppression of star formation along specific 
portions of spiral arms. Only Hdp exceeds the observed arm surface densities at radii 35"<i?<60", where 
gas depletion times are measured to be longest (Figure [1]) , as well as at larger radii (i.e. i?~80" and 
i?~ 100", where Td ep is also high.) As plotted in Figure [151 dynamical pressure is capable of supplying 
stability where other mechanisms fall short. 

Note that, even if with adjustments to the empirical calibration factors a and ola, cloud stabiliza- 
tion through Coriolis forces and shear cannot account for the differences in star formation properties 
between neighboring radial zones. Only dynamical pressure entails a non-monotonic radial variation in 
gas stability, given the pattern of streaming motions. 
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Fig. 14. — (Left) Comparison of the observed bar/arm (E{, a 



v j avmt 



solid black) and interarm (Ei 



dashed 



black) molecular gas surface densities in M51 with the critical surface densities for stabilization against gas self- 
gravity via coriolis forces (the so-called Toomre critical surface density Et oorare ; gray dotted), tidal forces (Etide', 
dashed gray) and shear [k replaced with Oort parameter A) in the disk (E s h ear ; red dotted) and in the arm 
(Esftear.sp; red solid). (Right) Comparison between the critical surface densities for shear, tidal and Coriolis 
forces (Toother', shown together as a gray band) with the surface density threshold for stabilization via dynamical 
pressure (Edp; dashed blue). The observed bar/arm surface density (E{, ar / arm ; solid black) is repeated from the 
left panel. 
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Fig. 15. — Critical surface densities vs. Eh 2 in the bar and spiral arms where different symbols and colors 
depict the potential for stabilization against gas self-gravity via Coriolis forces (the so-called Toomre critical 
surface density Et oomre ; gray open symbols), shear due to differential rotation (E 3 h ear ; open red) and including 
non-circular motions (E s ^ earjap ; solid red) and dynamical pressure (solid blue). Measurements are extracted in 
1.5" bins from the radial profiles in Figure [T?] The solid black line shows the line of stability. 
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